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ABSTRACT 


A  calculation  of  proton-proton  bremsstrahlung  is  presented 
that  is  fully  relativistic  and  gauge  invariant.  The  external- 
emission  part  of  the  amplitude  consists  of  a  sum  of  Born  one- 
boson-exchange  terms.  Elastic  cross-section  data  are  used  as 
input,  to  constrain  the  soft-photon  behaviour  in  accordance  with 
the  Burnett-Kroll  theorem.  The  problems  associated  with  the 
non-unitarized  one-boson-exchange  model  are  discussed  in  detail. 
Internal-emission  contributions  from  w  radiative  decay  and  A 
excitation  with  tt  and  p  exchange  are  added  coherently.  Computed 
cross  sections  are  compared  with  experiments  with  incident-beam 
energies  99  MeV  and  above.  The  u)  radiative-decay  contribution 
is  found  to  be  small  but  not  negligible  (<10%)  for  the  TRIUMF 
experiment,  and  may  have  to  be  considered  in  a  careful  comparison 
between  potential-model  calculations  and  the  experiment.  A  large 
discrepancy  with  the  asymmetric-geometry  data  from  the  Orsay 
experiment  is  noted,  similar  to  that  found  previously  in  potential- 
model  calculations.  The  A  excitation  contributions  are  included 
in  a  comparison  with  the  data  from  the  730  MeV  UCLA  experiment. 
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CHAPTER  1 


INTRODUCTION 

There  has  been  a  great  deal  of  work  done  in  the  field  of 
nucleon-nucleon  bremsstrahlung  (N+N N+N+ Y  )  since  the  paper  of 
Ashkin  and  Marshak,  published  in  1949  Z17-  Interest  in  this 
reaction  stems  mainly  from  the  possibility  of  learning  something 
about  the  off-shell  behaviour  of  the  N-N  force.  In  particular, 
it  is  hoped  that  the  reaction  can  serve  as  an  off-shell  testing 
ground  for  the  various  N-N  potentials  and  other  models  used  in 
the  treatment  of  N-N  scattering.  Other  reactions  involving  two 
nucleons  and  at  least  one  other  particle  (e.g.  pion  production) 
might  also  serve  this  purpose,  however  the  case  where  the  third 
particle  is  a  photon  seems  to  be  the  simplest,  since  it  avoids 
the  complications  of  three  strongly  interacting  particles. 

We  should  note  that  the  correct  off-shell  N-N  interaction  is 
important  to  our  understanding  of  the  properties  of  nuclear  matter 
and  of  other  many -body  systems  involving  nucleons. 

The  first  extensive  calculation  of  proton-proton  bremsstrah¬ 
lung  (henceforth  abbreviated  as  pptf)  was  done  by  Sobel  and 
Cromer  /  2_7  in  "1963-  Soon  thereafter,  the  first  experiments  were 
done  at  Harvard*  by  Gottschalk,  Shlaer,  and  Wang  £hJ  with  158  MeV 
incident  protons,  at  Manitoba  by  Warner  [_  5/  with  Tlab  =  48  MeV, 
and  at  Rochester  by  Rothe,  Koehler,  and  Thorndike  [IJ  with 

T,  ,  = 204  MeV  incident  protons, 
lab 


*Bremsstrahlung  from  collisions  of  protons  on  complex 
nuclei  was  first  observed  in  1952  by  Wilson  Z17,  who  attributed 
the  X  rays  to  n-p  bremsstrahlung.  Here  we  are  mainly  concerned 
with  bremsstrahlung  from  protons  on  targets  of  free  protons. 
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In  the  Harvard  and  Manitoba  experiments  the  directions  and 
energies  of  the  two  final-state  protons  were  measured,  while  the 
photon  was  undetected* *;  such  an  arrangement  is  now  commonly  re¬ 
ferred  to  as  Harvard  geometry**.  Since  the  A/ (\/ V  system  has  3 
degrees  of  freedom  -  9  momentum  components  minus  4  energy-momentum 
constraints  -  Harvard  geometry  once  overdetermines  each  event. 

In  the  Rochester  experiment,  on  the  other  hand,  the  directions  of 
all  three  final-state  particles  and  the  energies  of  the  protons 
were  measured,  thus  overdetermining  each  event  three  times. 
Nevertheless  the  experiments  are  difficult  because  of  large  back¬ 
grounds  and  low  counting  rates  -  the  ppV  cross  section  is  typ- 
-4 

ically  10  times  the  p-p  elastic  cross  section. 

There  have  been  several  experiments  performed  since  the 
pioneering  ones  mentioned  above.  Jovanovich  DJ  has  classified 
the  experiments  as  first  or  second  generation;  the  former  meas¬ 
ure  the  angles  of  the  protons  through  collimation  in  proton  tel¬ 
escopes,  while  the  latter  use  coordinate-measuring  proton  detec¬ 
tors  such  as  wire  chambers,  allowing  a  larger  region  of  phase 
space  to  be  simultaneously  observed.  An  important  problem  in 
the  analysis  of  the  data  is  the  correction  for  finite  solid-angle 
acceptance  (the  acceptance  window)  of  the  detectors.  This  is 
less  of  a  problem  in  the  second  generation  experiments  since  the 
acceptance  window  can  be  arbitrarily  set  in  the  computer  program 

V 

*The  Harvard  group  used  a  Cerenkov  counter,  but  only  to 
confirm  the  presence  of  the  bremsstrahlung  photons. 

**The  term  ’’Harvard  geometry”  is  used  in  a  more  restrictive 
sense  by  some  authors.  Our  usage  of  this  term  is  the  same  as 

that  of  Liou  DJ- 


when  the  data  are  analyzed. 


Besides  the  Rochester  experiment  there  have  been  three  other 
second  generation  experiments:  the  42  MeV  measurements  at  Manitob 
/  9_7  and  200  MeV  at  Triumf  /_  10_7i  both  in  the  Harvard  geometry, 
and  the  730  MeV  measurement  at  UCLA*  /_  1 1_/  in  the  Rochester  geome¬ 
try. 

The  data  from  Harvard  geometry  experiments  are  usually  pre¬ 
sented  in  the  form  of  distributions  dcr  / cl  sr  -  r)  vs 

(where  the  polar  axis  is  along  the  incident  beam  and  the  sub¬ 
scripts  3i  4,  Y  refer  to  the  scattered  protons  and  the  photon), 
and  integrated  cross  sections  ^ /dSl-^  (/tb/sr*)  .  The 

Rochester  group  presented  its  data  in  the  form  dcs/d Slcyy,  d J7r  d  k  , 
where  <^J7Cvy,  is  the  solid-angle  element  of  "p^  -  ~p^  and  k  is  the 
photon  energy  -  all  quantities  being  given  in  the  center-of- 
momentum  frame.  The  730  MeV  UCLA  data  are  reported  in  the  lab 
frame  as  spectra  d<r/d J23  d cl  k  in  k. 

For  a  more  detailed  discussion  of  the  experiments  performed 
so  far  the  reader  is  referred  to  two  review  articles:  one  by 
Halbert  £/\2.J  covers  the  period  up  to  1971  »  the  other  by 
Jovanovich  /~8  7  covers  the  period  since  then  till  1978.  Table 
1.1,  compiled  from  these  two  papers,  is  a  listing  of  the  ppY 
experiments  performed  up  to  now  with  beam  energy  99  MeV  and 
above • 


*In  this  experiment  the  directions  of  all  three  final-state 
particles  and  the  energy  of  one  of  the  protons  is  measured;  each 
event  is  thus  twice  overdetermined. 


Beam 

Energy 

(MeV) 

Institution 

Publication  Dates 
(first  -  last) 

Ref. 

99 

McGill  Univ. 

68  -  70 

48 

156 

Orsay 

72 

54 

158 

Harvard  Univ. 

65  -  67 

4 

200 

TRIUMF 

77 

10 

204 

Univ.  of  Rochester 

66  -  67 

6 

730 

UCLA 

77 

11 

Table  1.1  ppy  experiments  with  beam  energy  99  MeV  and 
above  (from  refs.  8  and  12). 

1.1  Potential  Model 

Most  of  the  calculations  have  been  done  in  the  framework  of  the 
potential  model.  It  is  hoped  that  differences  in  the  cross  sections 
predicted  with  various  N-N  potentials  will  be  larger  than  experimental 
errors,  thus  allowing  us  to  rule  out  some  of  these  potentials.  Unfor¬ 
tunately,  in  spite  of  the  number  of  experiments  and  calculations  this 
condition  has  never  been  sufficiently  satisfied.  The  reasons  are  a 
combination  of  the  following  factors:  1.  experiments  have  been  performed 
in  kinematic  regions  where  model  dependence  is  small;  2.  uncertainties 
due  to  finite  solid-angle  detector  acceptance;  3.  large  experimental 
errors;  4.  uncertainties  in  the  predicted  cross  sections  due  to  the 
neglect  of  some  corrections,  such  as  relativistic  effects,  coulomb 
effect,  rescattering  terms,  gauge  terms,  meson  currents,  etc.  Another 
problem  is  that  the  existing  N-N  potentials,  such  as  the  Hamada- 
Johnston  [  67  _7,  Reid  Z"68_7,  Bryan-Scott  ZT  69  _7 ,  etc.,  are  not 
really  equivalent  on  the  energy  shell.  These  on-shell  differences 
must  somehow  be  taken  into  account  before  purely  off-shell  differences 
can  be  extracted.  A  common  method  for  this  correction  (used  for 
example  in  ref.  13  and  14)  has  been  to  consider  the  ratio  of 


the  cross-section  computed  with  a  given  potential  to  the  cross- 
section  computed  with  the  same  potential  in  the  on-shell  approx¬ 
imation  (discussed  below). 


The  formulation  and  pioneering  calculations  of  pp*  using 
potentials  were  done  by  Sobel  and  Cromer  /2J7-  The  formalism 
starts  with  the  distorted-wave  Born-approximation  (DWBA)  ampli¬ 
tude 


i  *?'>  (i.D 

where  is  the  exact  outgoing  (incoming)  wave  function 

GUI 

with  the  nuclear  potential,  V^, alone*.  V  is  usually  taken  as 


/(«r»0  ,  — >,N  T? 

v  (p,o)  - 


+  L 


(t 


*t) 


(1.2) 


where  e,  jm,  and  p  is  the  charge,  magnetic  moment,  and  momentum 
of  the  proton,  and  and  k  is  the  polarization  and  momentum  of 
the  photon.  Gauge  terms  arising  from  the  replacement 
in  (whether  is  implicitly  momentum  dependent  as  from 
orbital  angular  momentum  terms,  or  explicitly  as  is  the  case  of 
the  Bryan-Scott  potential)  have  been  neglected.  Relativistic 
corrections  to  Vem  have  been  considered  by  Liou  et  al  /15,  l6/« 
The  amplitude  then,  is  given  to  first  order  in  the  electromag¬ 
netic  interaction  and  to  all  Born  orders  in  the  strong  inter¬ 
action. 


The  matrix  element  (1.1)  can  be  rewritten  as 

OO  \  ,  re*"') 


T* 


=  \/K‘~'> 
Vfn 


4* 


rc*  4-  . 


*  it-  &  V 

,  rM  <'»•>  re*>  , 

Y\  ^  r»  V  rr»  KV\ 


(1.3) 


m 


*The  Coulomb  potential  -  in  the  case  of  pp^  -  has  been 
neglected  in  most  calculations;  it  has  been  considered  by 
Signell  and  Marker  /l8_/  and  by  Heller  and  Rich  /  I?/. 
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where  =  (S  ~  Ha  -t  Z  is  the  free  propagator  and 

Ui  =  <<*>„  I  VM  I  /Tffw>  is  the  off-energy-shell  two-nucleon  matrix 
element.  The  first  term  in  (1.3)  represents  nuclear  scattering 
followed  by  photon  emission,  the  second  term  represents  photon 
emission  followed  by  nuclear  scattering,  and  the  third  term  is 
the  rescattering  (or  double  scattering)  term.  These  processes 
are  diagrammed  in  fig.  1.  The  rescattering  term  has  been  in¬ 
cluded  by  V.R.  Brown  /19_7»  Drechsel  and  Maximon  /20_7»  and  by 
Heller  and  Rich  CvJ- 

The  multitude  of  published  similar  potential-model  cal¬ 
culations  is  partly  a  result  of  some  serious  discrepancies  in  the 
early  calculations.  The  cross  section  predicted  by  Sobel  and 
Cromer  /"2_7  was  about  four  times  the  experimental  value  at  158 
MeV  measured  by  Gottshalk  et  al  /"4a _/,  and  ten  times  the  exper¬ 
imental  value  at  48  MeV  measured  by  Warner  C\J-  The  calculation 
was  done  in  the  lab  frame  (and  in  the  transverse  gauge,  •  k  ~  0  ) 
and  neglected  the  rescattering  term.  Another  calculation  was 
done  by  Duck  and  Pearce  based  on  the  Sobel  and  Cromer 

formulation  and  using  the  Tabakin  potential  Z?°_7-  Again  the 
rescattering  term  was  neglected  and  the  transverse  gauge  was 
used,  but  this  time  the  amplitudes  were  evaluated  in  the  CM 
frame.  The  results  agreed  with  experiment  at  158  MeV  but  were 
small  by  a  factor  of  three  at  50  MeV.  Subsequently  an  error 
was  discovered  in  the  Sobel  and  Cromer  equations.  Now  the  Duck 
and  Pearce  calculation  agreed  qualitatively  with  both  experiments, 
while  the  Sobel  and  Cromer  cross  sections  were  approximately 
correct  at  1 58  MeV  but  were  still  a  factor  of  6  too  large  at 
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48  MeV.  Also  there  were  problems  with  Sobel  and  Cromer's  dif¬ 
ferential  cross  sections  <j  9*  at  158  MeV;  at  48  MeV 

they  predicted  a  dipole  distribution  (although  Warner  did  not 
obtain  differential  cross  sections  from  his  experiment). 

This  discrepancy  was  resolved  by  Signell  /  22_/  who  derived 
an  analytic  expression  for  the  differential  cross  section, 

dOv  in  the  very  low  incident-energy  limit  (scattering- 
length  limit).  He  found  in  the  CM  frame  an  electric-quadrupole 
contribution  from  the  single  scattering  terms,  while  the  rescat¬ 
tering  term  was  negligible.  In  the  lab  frame,  on  the  other  hand, 
the  single  scattering  terms  have  a  large  electric-dipole  compo¬ 
nent  but  this  is  cancelled  by  a  large  rescattering  term.  This 
cancellation  takes  place  in  all  Born  orders  (except,  of  course, 
the  lowest,  where  the  rescattering  term  does  not  contribute). 

In  the  lab  frame  then,  the  rescattering  term  is  large  even  at 
low  energies,  and  must  be  included.  Subsequent  potential-model 
calculations  were  all  done  in  the  CM  frame.* 

1.2  Model-Independent  Calculations 

Another  method  of  obtaining  bremsst rahlung  cross  sections 
is  to  make  an  approximation  such  that  only  the  on-shell  nucleon- 
nucleon  scattering  amplitudes  are  used.  This  is  what  is  meant 
by  "model  independence":  although  the  on-shell  amplitudes  may 
be  obtained  from  a  model,  the  radiative  amplitude  is  to 


*It  is  interesting  to  note  that  before  this  error  (neglect 
of  the  rescattering  term  in  the  lab  frame.)  was  known,  an  early 
UCLA  experimental  group  of  Slaus  et  al  /  23./  measured  distri¬ 
butions  in  dy  which,  they  stated,  were  in  agreement  in  shape 
with  the  lab- frame  calculation;  i.e.,  their  distributions  appear 
to  have  the  incorrect  dipole  shape. 


» 


some  approximation,  uniquely  specified. 


Model-independent  cal¬ 


culations  are  useful  for  exploring  kinematic  regions  and  for 
checking  model  calculations.  The  approximations  that  have  been 
used  are  of  two  types:  the  two-energy  approximation  and  the 
soft-photon  approximation.  We  discuss  each  in  turn. 

Two-Energy  On-Shell  Approximation 

One  of  the  important  results  of  Sobel  and  Cromer's 
analysis  ^a4_7  iS  the  angular-momentum  expansion  of  the  off- 
shell  nuclear  matrix  element,  I  1  A^c+')/>,  in  terms  of 

the  phase  shifts  and  what  they  called  the  "quasiphase  parameters" 
This  is  entirely  analogous  to  the  phase-shift  expansion,  of 
Stapp,  Ypsilantis,  and  Metropolis,  for  the  elastic  scattering 
matrix  elements.  For  the  uncoupled  states*  the  expansions 
are 

=  2.  4  Al  (pn,  p/)  (1.4a) 

L. 

i  L  C  P-^) 

=  2.  (2L-f\)  Pu(u>s9)  -Z  AAp^Pf)  (1.4b) 

u 

Here  the  matrix  element  t  .  (t„  )  is  written  in  the  initial 

m  in 

(final)  nucleon-nucleon  CM  frame,  P^  Pf,  Pn  are  the  CM  nucleon 
momenta,  SL( p)  is  the  phase  shift, 

(p^J  Pc^  ~  ~yn  ^  r)  VtoCr)  (pc  jp)  rdr  0.5) 

o 

is  the  quasiphase,  and  UuCp.r)  is  the  radial  wavefunction  with 


*To  avoid  the  complications  of  nucleon  spins,  in  these 
introductory  remarks,  the  coupled  states  are  not  considered. 
The  equations  including  the  coupled  states  and  other  details 
are  given  in  the  paper  by  Sobel  ^  24_7 • 
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the  usual  asymptotic  condition: 

U.L  C  pj  r3  >  S i  ^  (  p  r~  -  ^  +-  Si.  C  p')') 

On  the  energy  shell  p  =  p.  in  t  .  and  p  =  p.  in  t>  ,  and  the 

n  *1  ni  fn  f  fn 

quasiphase  is  simply 

A/.  (  Pj  p)  =  S/’n  $l  Cp)  (1.6) 

Within  the  framework  of  the  potential  model  the  problem 
then  is  to  evaluate  the  integrals  (1.5)  for  the  quasiphases. 
Since  the  model  dependence  is  all  in  the  quasiphases,  it  can  be 
investigated  by  comparing  quasiphases  evaluated  with  different 
potentials,  although  any  model  dependence  found  in  this  way  must 
be  considered  as  preliminary,  since  the  total  matrix  element  may 
involve  cancellations.  Often  the  quasiphases  are  divided  by 
their  on-shell  values  (1.6)  before  the  comparison  is  made,  in 
order  to  remove  on-shell  differences.  The  results  are  that  at 
incident  lab  energies  less  than  about  100  MeV  the  model  depend¬ 
ence  is  small  and  not  likely  to  be  seen  in  the  present  experi¬ 
ments. 

The  two-energy  on-shell  approximation  consists  of  the 

replacement  of  the  quasiphases  by  their  on-shell  values  (the 

rescattering  term  being  ignored).  The  nuclear  matrix  elements, 

t  .  and  t„  ,  then  involve  only  the  phase  shifts  at  the  momenta 
ni  fn 

p^  and  p^  respectively.  This  approximation  has  been  studied  by 
Signell  and  Marker  /T8,  22,  257-  They  found  it  to  be  accurate 
at  small  incident  energies,  even  at  small  proton  exit  angles  - 
where  the  photon  energies  are  largest,  with  the  error  rising 
rapidly  with  increasing  incident  energies.  For  example,  at 
20  MeV  and  proton  exit  angles  9  =  35°  (coplanar,  symmetric 
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Harvard  geometry)  they  found  the  error  in  ^Ull)  dJlH  to  be  less 
than  1%,  rising  to  10%  at  61.7  MeV  9  =  30°  and  about  15%  at 
99  MeV  9  =  30°.  These  results  discourage  us  from  looking  for 
off-shell  effects  in  experiments  performed  at  low  energies. 

Finally  it  should  be  mentioned  that  the  two-energy  approx¬ 
imation  is  useful  for  examining  the  effect  of  the  on-shell  dif¬ 
ferences  among  the  potentials  and  for  correcting  for  these  dif¬ 
ferences  . 

Soft-Photon  Approximation 

The  second  type  of  model-independent  calculation,  the  soft- 
photon  approximation,  is  based  on  Low’s  theorem  /26_/.  For 
photons  of  low  energy  k,  the  radiative  amplitude  has  a  1/k 
singularity  (resulting  in  the  well  known  "infrared  catastrophe") 
and  can  be  expanded  in  powers  of  k  as 

T  =  ^  +-  T.  -  T,  k  *  0<V)  (1.7) 

Low’s  theorem  states  that  the  leading  two  terms  in  this  expansion, 
T  ^  and  T  ,  depend  only  on  the  non-radiat ive  amplitudes.  The 
first  term,  T  ^ ,  is  proportional  to  the  on-shell  amplitude,  while 
Tq  is  proportional  to  the  derivatives  of  this  amplitude  (with 
respect  to  the  two  invariants  of  the  non-radiat ive  process). 

The  essential  part  of  this  theorem,  that  Tq  is  model-independent 
as  well  as  T  ^  is  a  consequence  of  gauge  invariance.  The 
explicit  expressions  for  the  coefficients  T_1 ,  Tq  were  derived 
for  N-N  scattering,  by  Nyman  /  27_7*  Nyman’s  elastic  amplitudes 
are  given  at  the  Mandelstam  variables  5  -  [(p,  +pi)  +  (p3^ 

t  -  u?i-pS*  -P’il/2- 


and 


The  choice  of  the  on-shell 
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point  has  been  shown  by  Fearing  [_  28  J  to  affect  the  radiative 
amplitude  only  in  0(k)  or  higher*,  as  long  as  this  on-shell  point 
has  the  correct  limit  as  k  approaches  zero. 

In  discussing  bremsstrahlung  it  is  useful  to  distinguish 
between  the  external-radiation  processes  -  nucleon-pole  terms  with 
the  photon  being  emitted  from  an  external  nucleon  line  (the  single 
scattering  terms  of  fig.  1),  and  internal-radiation  processes  - 
all  the  others.  The  latter  type  are  represented  in  fig  2.,  and 
include  the  rescattering  terms  of  fig.  1,  meson-current  terms, 
and  isobar  excitation  terms. 

With  Nyman's  choice  of  the  on-shell  point,  the  leading  term 
in  the  expansion  (1.7)  is  solely  due  to  external  radiation  and 
is  gauge  invariant  by  itself,  the  1/k  singularity  being  due  to 
the  nucleon  propagator  in  the  radiating  leg.  The  second  term,  TQ, 
contains  a  part  from  the  external  processes,  which  by  itself  is  not 
gauge  invariant,  and  a  part  from  internal  processes.  This  internal 
part  is  independent  of  k^  and  makes  the  second  term  gauge  invariant 
as  well. 

In  the  soft-photon  approximation  the  two  leading  terms  in 
the  k  expansion  are  obtained  from  the  elastic-scattering  amp¬ 
litudes  while  higher-order  terms  are  irrelevant.  The  cal¬ 
culations  were  first  done  by  Nyman  [  27  J  who  used  a  numerical 
representation  for  the  amplitudes  to  obtain  the  cross  section. 

The  method  has  the  advantage  over  the  two-energy  approximation 
of  Signell  in  that  it  is  gauge  invariant  and  fully  relativistic, 

*More  precisely,  the  first  two  terms  in  the  expansion  (1.7)  may 
individually  change  with  the  on-shell  point,  but  their  sum  must 
change  only  in  0(k)  or  higher. 
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and  is  therefore  more  useful  for  the  higher  energy  experiments 
(Tfab  >  100  MeV)  which  are  of  most  interest  to  us. 

Squaring  (1,7)  gives  the  following  expansion  for  |H  : 


T-,  I*)  ,  , 

k  L|T„|  +2H<(T-,  l,  )+0Zle)(i.8) 


The  calculation  of  the  first  two  terms  is  greatly  simplified  by 
a  theorem  of  Burnett  and  Kroll  Z29_7-  They  showed  (for  fermion- 
boson  scattering)  that  for  unpolarized  cross  sections  (i.e. 
fermion  spins  unspecified)  the  leading  term  is  proportional  to 
the  non-radiat ive  cross  section  and  the  second  term  is  propor- 
tionalto  the  derivatives  of  the  non-radiat ive  cross  section. 

The  explicit  formula  for  N-N  scattering  was  derived  by  Fearing 
who  also  obtained  a  generalization  of  the  Burnett-Kroll 
theorem  for  polarized  cross  sections  (i.e.  nucleon  spins  being 
specified)/  28J7 • 

The  complete  model-independent  calculation  requires  the 
inclusion  of  the  |  Ti )  term  in  (1.8)  as  well  as  the  two  leading 
terms  discussed  above.  In  certain  kinematic  regions  neglect  of 
this  term  may  even  result  in  negative  "cross  sections".  The 
importance  of  the  various  model-independent  terms  depends  on  the 
incident  energy,  as  discussed  by  Fearing  (_  31_7*  We  shall  return 
to  this  later. 

1.3  Relativistic  Boson-Exchange  Model 

In  this  thesis  we  have  computed  the  pptf  cross  sections 
directly  from  relativistic  perturbation  theory;  the  strong 
interaction  was  generated  by  exchanges  of  mesons.  In  com- 
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parison,  the  potential  model  has  the  advantage  of  treating  the 
strong  interaction  to  all  orders  but  suffers  the  disadvantage  of 
being  non  relativistic  -  and  one  expects  sizable  relativistic 
corrections  even  at  158  MeV  incident  lab  energy.  Furthermore, 
none  of  the  potential -model  calculations  is  gauge  invariant, 
since,  as  mentioned  before,  gauge  terms  from  the  momentum- 
dependent  parts  of  the  potential  are  neglected. 

The  present  calculation  is  along  similar  lines  to  that  of 
Baier,  Kiihnelt,  and  Urban  Z  52.7.  It  has  the  advantage  over  the 
potential  method  of  being  fully  relativistic,  gauge  invariant  and 
automatically  includes  processes  such  as  that  in  fig.  3»  which 
contain  a  virtual  antibaryon  in  the  intermediate  state.  It  also 
suffers  from  a  serious  disadvantage:  the  strong  interaction  can 
only  be  handled  easily  in  lowest  order,  inclusion  of  processes 
of  higher  order  in  the  strong  interaction  is  impractical  with 
this  method.  Unfortunately  this  difficulty  does  not  appear  to 
have  been  fully  appreciated  by  Baier  et  al ,  with  the  result  that 
their  model  fails  to  conform  to  the  Low’s  theorem  calculation  in 
soft -photon  regions;  such  an  agreement  is  mandatory  for  all 
models.  This  difficulty,  and  the  method  we  used  to  handle  it, 
is  discussed  later  in  detail. 

We  have  also  studied  contributions  from  isobar  excitation 


and  meson  currents. 


CHAPTER  2 


KINEMATICS  AND  PHASE-SPACE  CONSIDERATIONS 


2.1  Kinematics 

We  shall  write  the  pptf  reaction  symbolically  as  1  +  2 -*3  +  4  +  f 
where  particles  1  and  2  are  the  beam  and  target  protons  respect¬ 
ively,  particles  3  and  4  are  the  scattered  protons  and  *  is  the 
emitted  photon.  For  convenience  the  photon  will  also  be  denoted 
as  particle  5-  Define  a  coordinate  system  in  the  lab  frame  with 
the  z  axis  along  the  beam  and  the  x-z  plane  in  the  horizontal 
plane.  In  this  system  the  incident  and  target  protons  have 
energies  and  momenta  t{  =  yr)  +•  ~p,  =  C  *  TU")]  '  2, 

~p3.-0  .  We  consider  only  events  in  which  the  two  final- 
state  protons  emerge  on  opposite  sides  of  the  y-z  plane;  the 
proton  with  a  positive  x  component  of  momentum  is  labelled  as 
particle  3,  that  with  a  negative  x  component  of  momentum  as 
particle  4.  In  terms  of  spherical  coordinates  we  have 


~pL  -  pc  (s  m<9t  c  a  +■  -Sin  9c  Sin  '/V  cj  ^CosBcz)  i  -  ^  %  5  (2.1) 

For  non-coplanar  events  the  momenta  are  sometimes  given  in 
terms  of  the  Harvard  angles  9C  ,  %  ,  defined  as  follows:  9c  is 
the  angle  between  the  projection  of  pc  onto  the  horizontal  plane 
and  the  beam  axis,  fj  is  the  angle  of  pc  out  of  the  horizontal 
plane.  Normally  6t  is  allowed  to  span  the  full  circle,  i.e. 

0  6  &v  £  2tt  ,  with  0<  0r  <  7T  for  kx  positive;  while  $3 , 
are  restricted  by  the  labeling  convention  to  half  circles, 
i.e.  0  —  @2>j  Qy-  ~  7T  .  (see  fig.  4). 
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In  terms  of  the  Harvard  angles  the  momenta  are 


PS 

p3  (  5 1  vo  $3  cos  ^3  3^  Si  to  V3  ^ 

4- 

£<C?5  $3  ^ 

ft 

^  (~"5 1  to  Ocj  CosV^y  "h  5m  Yc,  lj 

4 

C  OS  P'i  Cos  ft,  2r 

(2.2) 

T 

=  k  (  SfVt  Cos'f'g  -  S in  fy 

4- 

Cos  cos'Pf'  2-'') 

The  transformation  from  spherical  to  Harvard  angles  can  be 
obtained  by  equating  the  components  of  (2.1)  and  (2.2).  For 
k  we  have 


S  1  C  C>S  P)r 

5  1  VO  9y  Cosier 

(2.3x) 

S/h  'f’r 

~  ~  5i|o@y  Si  1 0  ^5" 

(2.3y) 

Cos  (9a-  cospy 

~  CoS 

(2.3z) 

therefore 

”  S.' vo  '  ("  “St  vo  #r  SmO  -5  — 

@7  =  polar  angle  of  the  vector 

Z  Z  M- 
'r  —  2. 

fa)  ,  ccC)  ”  f  cos  9r 

1  •  ©  •  y  9  X  h*  (\Y^ 

j  5  /  n  cos  J 

6  6  a  VO  &y  COi  4V") 

(2.4) 

taken  in  the  appropriate  branch  . 


Similarly 


^2,  5 1  y\  (  5 1  h  $3  5 1  h  fO 

02  ~  ~Lav\  1  (  cos  9^) 


-11  Z 


c  ft  £ 


77 


0  C  &2  4  - 


(2.5) 


and 


%  5 I  \o  (  S ih  ^  n'Pc'l 


TL  f.  Z  JL 

2.  —  -  2_ 


Ou  s  '  C  tan  Cosfii')  0  7=  6c,  \ 


(2.6) 
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As  a  measure  of  noncoplanarity  Gottschalk  et  al  /  4_/  de¬ 
fined  the  angle  <§  ~  ;  on  the  other  hand  Drechsel 

and  Maximon  /20_7»  who  used  spherical  coordinates  defined  the 
angle  -  ~k  L’fs  ^  (9T  ~  ftyV]  •  When  ©^  =  ©^  =  ©  these  two 
angles  can  be  simply  related  as  follows:  rotate  the  coordinate 
system  about  the  beam  axis  so  that  f 3  -  7T  ~  .  Then  =  %  -  7T 

and  §  -  —  %  and  from  (2.5)  we  have 


S 1  n 


§ 


S  1 


9 


(2.7) 


Since  the  quantities  involved  are  invariant  under  rotation  about 

the  beam  axis,  this  equation  holds  even  when  f3  7-  jj  ~  ^ 

For  a  fixed  BY  the  noncoplanarity,  <£>  ,  has  a  maximum  value 

allowed  by  kinematics,  (  ®r)  •  When  0Y  is  unrestricted  the 

noncoplanarity  has  a  maximum,  ,  at  a  certain  gamma-ray 

— ^ 

momentum  k0  with  angles  Oy~B 0j^-%  “  the  "limiting  gamma  ray". 

We  shall  now  consider  the  energy-momentum  equations  explic¬ 
itly.  In  the  lab  system  they  are 


pi  +  F-  *  **  =  pi 

(2.8) 

E~ Z  +  ^  '*7  +  ^  ~ ' 

(2.9) 

where 

=  Vp^  -- 

£"40+  ~  ^  ^ 

For  convenience,  in  what  follows  we  shall  consider  the  coordinate 
system  defined  so  that  ^3  =  71  -  %  =  $  .  The  components  of 

(2.8)  in  the  x,  y,  z  directions  respectively  are 


(p3s,'n93  -  si"  0^  cos^  +  k  St  n6r  oos^r  -  0  (2.10) 

(  p  3  5/Vi  @3  +  p  <y  sin  5/'h^  f  k  sin  0r  Sin  fr  =  0  (2.11  ) 
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P3  03  +  cos  +  kcosQ#  -  p, 


(2.12) 


Squaring  (2.7)  and  (2.8)  and  eliminating  fV  we  get 

(  p3  5/10  $3  -  ply  s/h  =  k  &r  ~  P  Pi  Pv  •Sf*  Pj  5 /  h  Ay  S/*7fi  (2.13) 

Eqs.  (2.9)  and  (2.12)  can  be  solved  for  C  9^  » 

the  result  (see  appendix  A)  is 


pH  =  < 


when  cos79r~ 

k  (2.14) 


^  P  ~  CoSZPyN) 

2  /}  Cc>3  6>i 

~Acos9l1  +  ccs  Pr  T  A  ^  ( co s' 0^  - co s\ 50~j 

\  Cc,s.*64,  - 

when  coi  2  Of  ^  COS  Ay 

where  P  =  (  4"iot  -  £"2P  c<?s0y  -  Pp,  -  pj  cos^  . 

Eqs.  (2.13),  (2.14)  and  (2.9)  can  now  be  combined  to  give  <£>  as 
a  function  of  p3  and  9*  .  We  have  maximized  this  function 
numerically  with  respect  to  p3  and  &v  to  obtain  and  at 

fixed  9r  with  respect  to  p3  to  obtain  <£  ( Q^)  .  Table  2.1 

shows  and  the  angles  0T  ,  fY  of  the  limiting  gamma  ray  for 

several  energies  and  angles  0^,  ©^. 

When  the  directions  of  the  final-state  protons  are  spec¬ 
ified,  the  locus  of  points  satisfying  the  kinematic  equations 
(2.8)  and  (2.9)  depends  on  one  other  parameter,  which  can  be 
taken  as  the  photon’s  polar  angle  .  Since  this  parameter 
is  periodic,  the  locii  are  closed  curves;  these  are  normally 
displayed  in  the  plane.  Such  plots  are  commonly  used 

by  experimentalists,  and  their  properties  are  well  known1 


For 


*See  for  example  the  review  paper  by  Halbert  /12/. 
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e3 

e4 

© 
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30 
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6.12 

72.4 

-90 

30 

30 

9.27 

63.4 

-90 

40 

40 

2.41 

60.0 

-90 

12 

12 

46.27 

63.O 

-90 

16.3 

16.3 

30.12 

62.7 

-90 

10 

18 

37-80 

62.6 

-84.7 

25 

25 

15-43 

61.4 

-90 

30 

30 

10.38 

60.3 

-90 

35 

35 

6.25 

58.7 

-90 

40 

40 

2.57 

56.3 

-90 

20 

20 

31.96 

52.6 

-90 

Table  2.1 

Maximum  noncoplanarity 

Wfor 

several 

energies,  Tlab 

,  and  angles  ©3,  ©4- 

®0 

and  fo  are  the 

spherical 

angles 

of 

the 

limiting  gamma 

ray.  Labelling 

is 

such 

that  10°  <  %  <  310 e  ,  -lo°  < 

f3  <  ^0 

O 

• 

1 

n  -  ^1)  —  • 

All  angles  are  in 

degrees . 
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example  the  locii  of  constant  k  are  the  lines  +  T,  =  T_  ,  -  k, 

3  4  lab 

hence  the  plots  are  convenient  for  locating  regions  with 

maximum  photon  energy.  In  appendix  B  we  show,  for  coplanar 
geometry,  that  the  closed  curves  are  approximately  elliptical 
when  p/m  Us  i' n  C&i,  +  <4  l  ;  the  ellipses  distort  even  at  low 

energies  when  ©^  +  ©^  is  made  small. 

2.2  Cross-Section  Formulas 

The  (invariant)  T  matrix  is  defined  by 


S$i  ~  +  L  (Ilf)  §  Yp3  +  p4  +  h  -  p,  -p,') 


X 


r  ^  -i 

V 

f  kvi 

U27T)34  J 

_  £~,  2  1  k )  _ 

(2.15) 


4c 


According  to  the  prescription  given  by  Bjorken  and  Drell  /  33 _/ 

the  differential  cross  section  is 

*  cCH), 


d  cr  - 


4c  i  ^3p3  ^  k 


,>u 


I'TT,  “ltl|  2  |  k  I  4  t-ty 

where  M  -  (  p3  +  pH  +  fc> 

>4  -  c  p,  + 


(2.16) 


and 


and  T2  are  the  velocities  of  the  beam  and  target  proton 
respectively.  Eq.  (2.16)  is  the  starting  point  for  deriving 
the  expressions  for  the  various  phase-space  factors  we  shall 
require . 

Rochester  Cross  Section 

From  (2.16)  we  get 


dcr/dJl%  dsiic/k  -  (2ttY  K  -VJe;  £i 


P3*  g(  p3  g(3p^ 


£; 


*7 


E  I  TP  S (£1  * £"i 


2  k 


H 


(2.17) 


Spin-* 
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where  )  H  I  7~ I  represents  the  sum  of  J  /f  T  over  final-state 

Spins  T'- 

proton  spins  and  photon  polarization  and  the  average  over 
initial-state  proton  spins.  The  second  delta  function  in  (2.17) 


determines 


P 4  as 


Pm  Pl  k  p3 
and  -  f  C  pc  ~~k  ~  psf  +  hi11] 

"  t  (fc  - ' +  Ps  ~  2  P3  I'pc 

where  u  is  the  angle  between  ~p^  and  p\ 
V/e  now  need  to  evaluate 


"  k  I  Coso( 

-  k. 


Wo 


=] 


with 


$  (  Pi'))  d  p3 


■f  ( P3 )  ~  ^3  f  ^  H  k  ~  ^  i 


(2.18) 


o£ 


(2.19) 


0£ 


With  =  (,p32+wVLy and  given  by  (2.18)  we  get 
/ d  p3  -  Pj/t“3  +  (Pv  ~  /  Pc  ~  k  I 

and  ^  6  (T  -P  C p3^y>  p>3  -  ^  S  f-fO  d-P  = 


cA£ 
d  Pi 


a?3  l-fCp3^0 


p3  ^  3)  P  M 


-  “P3 


(2.20) 


It  is  understood  that  p^  is  a  zero  of  (2.19) • 


Putting  these  results  back  into  (2.17)  we  obtain 


&<S 


wv\ 


P/  H 


ft 


MidJi-tiV  v,|r,£-t  |  ( c5 * e,-) - ft) e; 

The  quantity  in  the  denominator  in  (2.21)  can  be  rewritten 
simply  in  terms  of  ~p^  as  follows: 


(2.21) 


i 
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I  Pi"  ( tj  *  fcn')  -  p3  ■  (f.  -T)r3 1  =  |  Pj  (h+h)  -j?3 •  | 

=  I  P3  ^  -  ft  •  ft  ft  | 


(2.22) 


In  the  lab  frame 


=  0 

Pt  =  ^ 
£"»  =■ 
or,  = 


Pi 


and  (2.21)  with  (2.22)  reads 


(2.23) 


d<j 


*>2i 


Jjli  clJly  dkl,ak  2.CStt)  p,  j  p/  t,  -p3.p*,  £TS 

In  the  center-of-momentum  frame 
'Pi  ~  -T^z  ~  ?,A~, 

pc  =  0 

and  (2.21)  reads 


j_ 
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(2.24) 


/  (2.23) 
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ol 42 -j  c4J2s  cl  k 
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k  pi" 


^2 

*  3p»hS 


^  ^Tp)  p,  t,  |  p3a  (h3  +  tO  +  pj  -T  fc: 


(2.26) 


In  the  experiment  of  Rothe,  Koehler,  and  Thorndike  /6  /  the 
differential  cross  sections  were  given  not  in  terms  of  ©^, 

T3,  but  9^,  ,  the  spherical  angles  of  the  vector  ~q  =  p^  -  "p^ 
(in  the  CM  frame).  In  this  case  the  cross  section  is  [2 7b/. 
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(2.2 7) 


AJl^dJlrdkj^  Pteri  p,  |^(t3  +  E-0-J2T-^ 

- 

Using  the  coordinates  of  q  instead  of  p  has  the  advantage  that 
in  the  soft-photon  limit  ©  becomes  the  usual  CM  elastic- 

q 

scattering  angle 
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Note  that  for  soft  photons  the  phase-space  factors  in 
(2.24),  (2.26),  (2.27)  can  be  expanded  as  a  power  series  in  k 
with  the  leading  term  of  0(k). 

Harvard  Cross  Section 

From  (2.15)  we  have  in  the  lab  frame 


<dc r _ 

ckJli  &Jb-\  cl  Qy 


yy\3  sm  9r  (  $C  Ypf  ~  pt')  !  "TJpJ  Pi  dp3  k  dkdfr  (2.28) 

zCz-iif  p,  £-3 


This  phase-space  integral  is  evaluated  by  Drechsel  and  Maximon 
/20_7;  the  result  is 


dp 


-J^3- _  p/~  »  y  I  J|z 

ET3  E1+  I  id  \ 


(2.29) 


where 


M  —  (p^Slyi0^  —  p  3  S  ir\  C  os  ( S  i  h  (Pj  +  ficj)  1X3  Slh  ZJiy  S 1  h 

~  k  (nr3  sir \  -TTij  S/'h@3  )  sird'dx 

+  2.  S/'n  63  S/rd£  [  p3  CO^  @3  ~  p<y  Coj  ~  (7^  p3-'Zrty  cos$y] 


It  is  understood  here  that  the  coordinate  system  is  defined  so 
that  %  =  77  -  f4  -  J  .  The  phase-space  factor  in  (2.29)  is 
singular  when  f .  This  is  easy  to  see  for 

symmetric  events,  =  9:  ^  fP^when  =  X  TT  » 

implying  (from  the  symmetry)  that  p^  =  p^  and  hence  that  N  van¬ 
ishes.  Furthermore  when  9v^0ori\  only  coplanar  events  are 
allowed,  hence  CO)  -  (t t)  =  0.  For  events  with 

0<  4>  <  (f^there  are  no  solutions  to  the  kinematics  when  0y  is 
such  that  •  For  example  consider  the  situation 

T  =  200  MeV,  ©,  =  6.  =  16.3°,  £  =  5°.  Fig.  5  is  a  plot  of 
lab  7 
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and  shows  that  tkvr,*?  f&v)~  £  at  &Y  ~  6°  and  at  9r-l6 3°  , 
hence  there  are  no  solutions  for  6V<6°and  for  Qr  >  163°.  Be¬ 
cause  of  these  two  problems  -  singularily  of  phase-space  factor 
at  the  kinematic  limits  and  restricted  range  of  Qr  -  Gottschalk 
et  al  defined  a  different  angle,  'Y V  ,  that  avoids  the  problems 
and  reduces  to  9r  in  the  coplanar  geometry.  Yr  is  defined  in 
terms  of  the  limiting  gamma  ray 

k0  -  leoSlV»0o  +•  kDcosO 0Z  (2.30) 

in  the  following  way:  consider  the  vector 

~k'  =  ~k  ~  U  ko  (2.31) 

where  d  is  a  scalar  such  that  has  no  y  component  (this  can 
always  be  done  since  ""£>  has  a  non-zero  y  component),  then  Yr  is 
the  angle  between  k '  and  the  z  axis;  i.e., 

=  k'  (2.32) 

The  angle  %  as  defined  by  (2.32)  can  range  from  0  to  2tt  .  (see 
fig.  6).  The  components  of  (2. 31)  read 

k  S'V 1  Yy  -  h  SioOy  Cos  fr  ~  o^s'fo  (2.33*) 


0  k  Sj  ^  Oj-  S/nYy  o(  ko  *0  0o  5/ to  Y0 


(2.33y) 


(2.33z) 


k'  COS  Yy  -  fe  ko  ODS0O 

From  (2.33y)  we  have 

o<  k0  =  Ck  5 /to  9r  S\\r\  /  (Y/n  £0  S)h  fo) 

Substituting  this  into  (2.33x),  (2.33z)  we  solve  for  Y"r  S 
=  k  Sih^r  -  -  5 infr  cot  f o') 

~  kCcos9r  ~  si\n9r  si»fr  tet&o  esc  ft) 

Yv  is  the  polar  angle  in  the  x1  -  x2  plane  of  the  point 
(Yt,  cO  -  f  cos  0*  -  s/m@Y  s/v>fY  cot0o  CSC  'Po,  S/’y)  9r('cos?r -Sih%  co-tfo)) 


(2.34) 
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The  phase-space  integral  of  the  Harvard  cross  section, 

d  ,  is  evaluated  by  Liou  and  Sobel  ^16_7  in 

terms  of  the  Harvard  coordinates  ©.  .  ?_•  defined  in  section  2.1. 

1  ’  1 

The  result  is  rather  lengthy  and  can  be  found  in  the  appendix 
of  their  paper;  in  terms  of  their  factor  F,  the  cross  section 
is  given  by 


d  cr 
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F  ^iTl2 

Sp/kyS 


(2.35) 


Elastic  cross  section 

We  shall  also  need  the  differential  cross  section  for 
elastic  N-N  scattering.  Here 
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(2.36) 


The  notation  is  the  same  as  before;  now  k  =  0. 
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In  the  CM  system 
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Integrating  over  the  delta  function: 

^  S(2f'-3f)p'Vp'  =  [  h{Z  Jp'%^2-  -2 e)  p'^p' 

Z  p /g!  £  <V  ~  p')  p'^dp' 


T  P 


UJ/ 


;  e'  —  e 


Putting  this  result  into  (2.38)  and  summing  over  the  spins  we 


get  the  expression 

d<r  ~  i  J 


dvPcm  C^TT)2-  P  £~Z 


(PI 
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5  pi** 


(2.39) 


c4 


Sometimes  the  invariant  d(-t)  is  required,  where  t  =  (p^  -  p^)' 
is  the  invariant  momentum  transfer.  In  the  CM  system 


i  ~  ~Lt  pX5i^x-~  -  ~  2  p1  ( I  - 


(2.40) 


where  9^  is  the  angle  between-^  and  "p*.  Since  cU2c*  ~  d0cvndfc 


C*-n 


and  because  of  symmetry  with  respect  to  the  azimuthal  angle 
dcr 


ci  Q  z.  *v» 

therefore 
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Substituting  from  (2.39)  we  get 
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CHAPTER  3 


ONE-BOSON-EXCHANGE  MODEL  CALCULATION  — 
EXTERNAL  EMISSION 


3>1  Introduction 

In  view  of  the  good  agreement  with  experiments  reported 
by  Baier,  Kiihnelt,  and  Urban  /"32 _7  (BKU),  we  decided  to  use 
their  method  to  obtain  a  basic  calculation  of  pp*.  We  could 
then  add  other  contributions  such  as  those  due  to  meson  currents 
and  A  excitation  as  refinements  to  their  basic  model. 

Although  we  were  able  to  reproduce  the  results  of  BKU  we 
found  their  model  to  be  in  gross  disagreement  with  the  data  of 
Nefkens  et  al  and  with  other  calculations  at  very  small 

photon  energies  (for  e.g.  at  T^^  =  200  MeV  9^=  9^=40°):  it 
predicted  cross  sections  that  were  much  too  large.  The  problem 
was  clarified  when  we  applied  the  same  model  to  elastic  N-N 
scattering:  the  S-partial  wave  was  found  to  be  in  serious  error 

due  to  the  neglect  of  unitarizat ion. 

In  our  calculation  we  have  attempted  to  correct  for  this 
deficiency  by  imposing  agreement  with  the  soft-photon  theorem 
as  a  constraint  on  the  model.  For  the  N-N  interaction  part  of 
the  amplitude  we  use  the  model  of  Erkelenz,  Holinde,  Machleidt 
(henceforth  cited  as  EHM)  /3 kj  in  the  Born  approximation  with 


*It  is  true  that  this  experiment  is  outside  the  range  of 
the  model,  because  of  the  large  730  MeV  beam  energy.  However 
the  calculated  cross  section  was  about  16  times  the  measured 
cross  section  even  at  small  photon  energies,  and  this  problem 
persisted,  though  not  at  the  same  magnitude,  at  lower  incident 
energies . 
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the  parameters  modified  to  satisfy  the  constraint  mentioned 
above . 

3.2  The  Model 

The  OBE  amplitude  is  a  sum  of  external-emission  diagrams 
with  one  nucleon  leg  off  mass  shell  (fig.  7)»  The  strong 
interaction  is  represented  by  the  lowest  order  exchange  of 
pseudoscalar,  scalar,  and  vector  mesons  with  form  factors 
included  at  the  meson-nucleon  vertices.  We  have  used  the  EHM 
OBE  model  as  a  guide  to  meson  parameters;  appendix  C  is  a 
summary  of  this  model.  For  more  detailed  information  the 
reader  is  referred  to  refs.  /3^,  35*  377  and  the  review  article 
by  Erkelenz  £j>&7 • 

In  our  pp/  calculation  we  have  tried  to  use  the  mesons, 
couplings,  form  factors  and  coupling  constants  as  given  by  EHM 
(table  C.1  and  eq.  C.3).  However  since  we  did  not  iterate 
our  N-N  T  matrix  through  the  Lippmann-Schwinger  equation,  as 
EHM  did,  modifications  had  to  be  made  to  restore  the  fit  to 
the  N-N  data.  We  used  the  criterion  that  the  pptf  cross  section 
must  have  the  correct  soft -photon  behaviour.  A  necessary  con¬ 
dition  then,  according  to  the  Burnett-Kroll  theorem  (section 
1.2),  is  that  the  model  predict  the  correct  elastic  cross 
section  and  its  derivatives.  This  condition  was  satisfied  by 
modifying  the  form-factor  parameters  (the  cutoffs)  and  by 
multiplying  the  ppif  cross  sections  by  a  correction  factor^ 

R  ,  equal  to  the  ratio  of  the  experimental  elastic  cross 
c 

section  to  the  computed  elastic  cross  section  at  the  on-shell 


point  (s,  t)  with 


(3.1) 
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We  found  that  as  long  as  100  MeV  <T_  ,  <  730  MeV  and  9~.,  was  not 

lab  CM 

too  close  to  the  Coulomb  peak,  the  computed  elastic  cross  sec¬ 
tions  had  approximately  the  correct  dependence  on  9^.  The 
correction  factor,  therefore,  can  only  have  a  weak  dependence 
on  9_.,  and  it  is  sufficient  to  evaluate  R  at  a  single  9^,. 

We  tabulated  R^  at  9^  =  60°  at  several  energies  from  9.&9  MeV 
to  730.5  MeV;  the  correction  factor  at  any  energy  was  then 
obtained  through  interpolation.  The  result  is  plotted  in 
fig.  8.  The  slight  waves  in  this  graph  are  artifacts  of  the 
cubic-spline  interpolation*.  The  cross  sections  were  computed 
with  the  T  matrix  given  by  eqs.  (c.4)  and  (c.5).  The  param¬ 
eters  used  are  those  of  table  C.1  with  the  exception  of  cut¬ 


offs,  which  are 


A-n  -  -  As  -  1 2  50  MeV 

Aa  -  100  0  MeV 
A  p  ~  Alo  ■  A$  6  50  MeV 
A*=  P-500  AleV 


}  (3.3) 


These  are  just  the  cutoffs  taken  from  EHM  (table  C.l)  scaled 


*Since  the  points  contain  statistical  uncertainties  curve 
fitting  would  be  more  appropriate.  Interpolation  was  used 
because  it  is  much  easier  and  gives  adequate  results. 
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by  a  factor  of  0.5  (with  the  exception  of  A,).  The  scaling 

factor  was  treated  as  a  parameter,  determined  simply  by  trial 

and  error  to  give  the  correct  dependence  of  d6/dflc„  on  9™,  and 

CM 

T_  V* 
lab 

The  Coulomb  interaction  was  treated  approximately  and 
incompletely  by  adding  a  one-photon-exchange  diagram  with  the 
"dipole”  form  factors.  The  tfpp  vertex 

->  i  a Fx(^)  (3.4) 

where  the  F1  and  F 2  are  defined  in  terms  of  the  electric  and 
magnetic  form  factors 

(i)  -  F,  (  i)  +-  -5—  P*.  ( f) 

and  Gn(i)=  F,  (t)  *■  2m  Fl(£)  (3-5) 

with  Ge  and  G^  accurately  given  by  k6_/ 

(t)  _  /  t- _ \'Z 

Ge  d)  ~  |  h-  X  "  (  '  ~  O.Jl  6-eV*  )  (3.6) 

where  1  +  X  —  2. 79  is  the  proton  magnetic  moment  in  units  of 
nuclear  magneton. 

Fig.  8  shows  that  for  100  MeV  <T_  ,<730  MeV  and  ©  =  60° 

JL3.D  Url 

the  model  gives  a  reasonably  good  dependence  on  the  lab  energy. 
The  distributions  in  ©^  at  a  few  energies  are  shown  in  fig.  9» 
For  N-N  bremsstrahlung  the  invariant  T  matrix  was  defined 
by  (2.15).  We  factor  out  the  polarization  vector,  of  the  emitted 

p 

photon  to  obtain  a  four  vector  T  ;  i.e., 

T  =  T'"  (3.7) 

(For  notational  simplicity  we  shall  often  not  distinguish  between 
a  matrix  and  the  matrix  element).  We  use  the  meson-nucleon 
couplings  as  given  in  eqs.  (C.3).  For  the  sum  of  the  two 
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diagrams  of  fig.  10 

T P2;  1  P  =  -e  3^  Fw  (  (p^-pi'f) 

*ZjTM(p,.k)  +HV'>SCpl-Wp,.i)\u,  <3-8> 

"  M\'  tta 

where  the  index  ct  specifies  the  meson  exchanged;  the  index 
i  =  1,  ...,  4  in  the  spinors  specifies  both  the  spin  and  mo¬ 
mentum,  thus  u.  =  u(s . ,  p.); 

l  i  •*!  ’ 


=  (? 


Kvi 


(3.9) 


p  J  -  '  '  V  (>  pZ  - 

is  the  nucleon  propagator;  ('p')  is  the  photon-emission 

k  1  0*,  /O  kJ  C01  J 

vertex;  r\y  and  are  the  meson-nucleon  vertices  at 

points  A  and  B  of  fig.  10  respectively;  Z_A  Co^  is  the  meson 
propagator;  and  F<*  ( <&)  are  the  strong  interaction  form  factors 

defined  in  table  C.l. 

For  reasons  given  in  section  3*3 »  in  writing  down  the 
photon-emission  vertex  we  neglect  the  fact  that  one  of  the 

nucleon  legs  is  off-mass  shell.  In  this  static  limit 


l  cr 


(3.io) 


where  e%/^  is  the  anomalous  part  of  the  proton  magnetic 
moment  C**— 1#79)« 

The  meson  propagators  and  meson-proton  vertices  are 
given  below,  in  table  3*1  ♦  with  q'1*  =  (p^  - 
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Meson 

A 

A  c-) 

scalar 

1 

1 

1 

o  *  _  *• 

^  Kn  cs) 

pseudo- 

i-  Xs 

i  rs 

) 

scalar 

*  ^ 

%  ~  **  Cps) 

vector 

^  ^  ^ 
3-m 

~  8VJ 

Table  Definition  of  the  meson-proton 

vertices  and  the  meson  propagators 
used  in  eq.  (3-8). 


The  gauge  term  ^  in  the  vector-meson  propagator  does 

not  contribute  because  it  couples  to  a  conserved  current.  For 
example 

l*-Lf  U.X  ^  ~  be  (  p>cj  -  p 2 ,'S)  U.x  -  0 

since  both  spinors  satisfy  the  Dirac  equation. 

The  last  two  diagrams  shown  in  fig.  7  can  be  obtained 
from  (3.8)  by  the  transformation  s1P1^>s2P2  ,  *  while 

the  exchange  terms  required  by  the  Pauli  principle  can  be 
obtained  from  the  four  diagrams  considered  above  by  the  trans¬ 
formation  s^p^s^p^  and  a  change  of  sign.  The  total  T  matrix 

is  then 

T  =  6^Tm  =  2;  3,4) 

+  7^(8,  Ij  Hj 3")  ~  -  T£&,l;  3,v)J 


(3.11) 
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Because  of  the  large  number  of  terms  in  the  amplitude 
(3«31),  the  usual  trace  method  of  summing  1  T I  over  the  proton 
spins  is  not  practical.  We  evaluated  the  matrix  element  for 
each  of  the  16  proton  spin  states  (and  for  each  ja,  )  numerically 
and  simply  summed  the  mod-squares  of  these  complex  numbers. 

The  sum  over  the  photon  polarization  was  done  by  means  of  the 
relation 

=  -TMT*  (3.12) 

To  make  use  of  this  formula  it  is  necessary  that  the  amplitude 
be  gauge  invariant.  Then  1^-0  and 

(k°T°)(k°T°f  =c?-f ±  1ZT-T+ 

Since  kM  is  a  null  vector  this  means  that 

-y°  ~T  •  'j* 

i.e.,  tM  is  a  space-like  or  null  vector  and  the  right-hand 
side  of  (3.12)  will  always  be  a  non-negative. 

It  is  easy  to  show  that  each  of  the  four  terms  in  (3«11) 
is  individually  gauge  invariant.  For  i  =  3 j  =  1*2 


r  [  rM(Pi  -  m"*  5^-00^]^ 


—  Cc 


j  [/  (  Pi  +  X  *V»)  +  (ft  ?)  /]  U 

2  pc  •  k  -2p  r  k 
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((~/i  +  -  zpi-k )  +  + to) ■  fe) 

2  Pl  •  k  “  <2  Pj  ’  k 
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£  C-W"''0  +  AP'P/uj  =0  gee/ 


In  the  above  we  used  the  commutation  relation  for  the  gamma 
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matrices  ^  #  +  $  £  =  -2a. b  and  the  fact  that  the  spinors  satisfy 
the  Dirac  equation  (jf.  -  m)  u.  =  u.  (if.  -  m)  =  0. 

Note  that  since  is  an  invariant  and  we  do  not  work 

in  any  particular  gauge  it  makes  no  difference  whether  the 
matrix  element  is  evaluated  in  the  lab  or  CM  frame  (in  contrast 
to  the  potential  model  calculations) .  Also  the  gauge  condition 
provides  a  handy  numerical  check  -  in  particular  at 
0r -  0  ,  since  then  kM  =  (k,  0,  0,  k)  and  T°  =  T^. 

3.5  Discussion 

We  shall  first  discuss  the  problem  of  using  the  off  shell 
N-N  amplitude  in  Born  approximation.  As  mentioned  before,  EHM 
unitarized  the  meson-exchange  pole  terms  through  the  Lippmann- 
Schwinger  equation.  Since  we  avoided  this  procedure,  we  had 
to  apply  corrections  in  order  to  restore  agreement  with  the 
soft-photon  theorem. 

To  study  this  problem  in  more  detail  in  terms  of  partial 
waves  we  used  the  simple  four-pole  N-N  model  of  Arndt,  Bryan, 
and  MacGregor  (ABM)  /"38_7*  The  four  mesons  considered  were 
the  tt  ,  p  ,  w,  and  a  scalar,  T=0  meson  of  mass  of  about  400  MeV. 
The  L>1  partial-wave  amplitudes  were  then  unitarized  geo¬ 
metrically;  this  amounts,  for  the  uncoupled  states,  to  pro¬ 
jecting  the  Born  amplitudes  onto  the  unitary  circle.  Two 
methods  of  doing  this  were  used  by  the  authors:  1.  taking 
the  Born  amplitude,  BL,  as  the  real  part  of  the  unitarized 
amplitude  -  the  imaginary  part  being  determined  by  elastic 
unitarity,  i.e*  =  sin  SL  to s  &i_  »  2.  taking  the  Born 

amplitude  as  the  phase  shift,  i.e.  B^  —  8L  .  We  should  note 


- 
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that  this  geometric  unitarizat ion  is  only  valid  when  the  model 
amplitude  is  small  for  uncoupled  waves). 

In  this  (ABM)  model  the  coupling  constants  were  obtained 

from  a  fit  to  N-N  data  in  the  energy  range  from  25  to  350  MeV 

by  means  of  the  reduced  second-derivative  matrix  of  Arndt  and 

13  3 

MacGregor  /~39_7*  The  Sq,  and  ,  phase  shifts  were 

excluded  from  the  fit  and  were  treated  as  free  parameters. 

For  pp  scattering  then,  the  model  gives  the  L>  1  phase  shifts; 
the  S  wave  can  be  put  in  from  phase-shift  analysis  in  order  to 
compute  observables*. 

The  authors  found  the  unitarity  corrections  to  be  un¬ 
important.  This  result  was  used  by  Baier,  Kuhnelt,  and  Urban 
to  justify  the  use  of  the  Born  amplitudes  in  their  NN*cal- 
culation.  The  fallacy  here  is  that  the  S  wave  was  excluded 
from  consideration  by  Arndt,  Bryan  and  MacGregor  and  must  be 
supplied  externally,  otherwise  the  result  can  be  a  catastrophic 
disagreement  with  data,  as  we  illustrate  below. 

d  <r_ 

Fig.  11  shows  the  elastic  cross  sections,  dJlc™  , 
evaluated  at  three  energies.  The  calculations  were  done  by 
projecting  out  the  partial  waves  and  summing  these  according  to 
the  formulas  in  Stapp,  Ypsilantis  and  Metropolis  (SYM)  /~4o_/ ; 
the  series  was  continued  until  the  desired  accuracy  was  reached. 
The  expressions  for  the  QBE  partial-wave  amplitudes  are  quite 
lengthy,  but  have  been  written  out  in  a  number  of  papers  (see 


♦Although  to  minimize  the  chi-squared  one  should  use  the 
S  wave  obtained  from  the  reduced  second-derivative  fit. 
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for  e.g.  refs.  41*  and  42).  We  took  the  Coulomb  scattering 
into  account  by  adding  the  non-relativist ic  Coulomb  amplitude 
(as  given  in  SYM);  i.e.  we  set  the  Coulomb  phases  in  the  nuc¬ 
lear  part  of  the  total  amplitude  to  zero.  This  amounts  to 
neglecting  the  Coulomb  potential  in  the  region  of  strong  inter¬ 
action  and  is  an  acceptable  approximation  here,  because  we  are 
mainly  interested  in  the  region  >  90  MeV  and  large 

enough  that  Coulomb  scattering  is  small  and  we  are  looking  at 
effects  that  are  considerably  larger.  The  solid  curve  in 
fig.  11  is  the  result  with  no  unitarizat ion  and  with  the  S  wave 
taken  from  the  OBE  model  unmodified.  The  dashed  curve  is  the 
result  with  the  L2:  1  waves  geometrically  unitarized  according 
to  Bt  =  sin&L  cos£l  (see  appendix  E  for  more  details).  Both 
of  these  calculations  are  in  gross  disagreement  with  data  /f?1_7- 
The  effect  of  the  unitarizat ion  is  relatively  small  and  increases 
with  increasing  energy  and  decreasing  scattering  angle.  As  the 
energy  increases,  generally  the  L2.1  phase  shifts  increase  and 
unitarity  corrections  become  increasingly  important.  On  the 
other  hand  the  largest  partial  waves  -  the  P  waves  -  are  odd 
functions  of  cos  ©CM  and  therefore  make  a  vanishing  contribution 
as  ©CM  approaches  90°. 

Fig.  12  shows  differential  cross  sections  computed  at  the 
same  three  energies,  but  this  time  the  S  wave  was  put  in  from 

_  l  £>  (  r  f\  <-  N 

phase-shift  analysis  (according  to  o(0  -  -z  Si*  o(  o,)  ). 

*Eqs.  (A.7c)  and  (A.7d)  in  this  paper  contain  a  typographical 
error:  the  expressions  should  not  be  multiplied  by  2J  +  1  (as 

they  are),  rather  they  should  be  divided  by  this  factor. 


The  solid  curve  is  the  result  with  L_>  1  unitarized  as  above, 
while  the  dashed  curve  is  the  result  with  no  unitarizat ion 
(except  the  S  wave  which  was  unitary).  This  time  the  agree¬ 
ment  with  experiments  is  much  better  and  the  unitarizat ion 
makes  a  small  improvement  at  the  higher  energies.  The  remain¬ 
ing,  relatively  small  discrepancy  can  be  attributed  to  the 
approximation  in  the  handling  of  the  Coulomb  scattering  and 
errors  in  the  input  of  the  S  wave. 

The  above  study  demonstrates  that  the  problem  lies  in  the 
S  wave,  which  the  unmodified  meson-pole  model  gives  incorrectly* 
This,  of  course,  is  not  surprising  since  the  S  waves  are  known 
to  be  outside  the  scope  of  the  meson-pole  model.  We  should 
also  point  out  that  the  problem  gets  worse  with  increasing 
energy,  for  unlike  the  true  S  wave  amplitude,  which  goes 
through  zero  at  Tlab—  250  MeV,  the  model  S  wave  increases  with 
energy  and  continues  to  dominate  the  total  amplitude.  In  the 
OBE  model  of  EHM  the  S  wave  is  corrected  partly  by  adjusting 
the  cutoffs  in  their  form  factors,  but  mainly  by  unitarizing 
through  the  Lippmann-Schwinger  equation.  Thus,  unitarity 
corrections  are  small  when  the  S  wave  is  excluded  -  as  in  the 
ABM  model,  but  are  very  important  when  applied  to  the  total 
amplitude  -  as  in  the  EHM  model. 

We  shall  now  make  a  few  remarks  about  the  previous  OBE 
NNy'  calculation,  of  Baier,  Kiihnelt  and  Urban  /  32/ ,  since 
the  authors  appear  not  to  have  fully  appreciated  the  problem 


*This  was  confirmed  by  examining  the  partial-wave 
amplitudes  and  the  phase  shifts. 
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discussed  above.  Their  statement,  that  the  S-wave  phase  shifts 
(of  the  ABM  model)  are  within  two  or  three  standard  deviations 
of  the  tabulated  values,  is  incorrect,  because  the  S-wave  phase 
shifts  that  ABM  refer  to  are  those  from  the  reduced  second- 
derivative  matrix  Z~76_7;  i.e.,  it  is  not  a  model  phase  shift,  rather 

it  is  that  S-wave  phase  shift  which  along  with  the  model  phase 
shifts  gives  the  best  fit  to  the  data. 

It  is  not  surprising,  therefore,  that  Baier  et  al  found 
gross  disagreement  when  using  a  number  of  N-N  OBE  models  coupling 
constants.  When  used  in  Born  approximation  these  models  should 
give  the  correct  L>1  partial  waves  but  unpredictable  S  waves. 
Baier  et  al  used  the  coupling  constants  from  the  model  of  ABM, 
but  made  small  corrections  to  these  to  improve  the  agreement  of 
their  photon  angular  distributions  with  experiment.  However, 
we  found  the  elastic  cross  sections  and  the  partial-wave 
amplitudes  to  be  quite  sensitive  to  these  variations  in  the 
coupling  constants.  Fig.  13  shows  the  effect  of  the  changes 
in  the  partial- wave  amplitudes.  The  dashed-dotted  curves  in 
fig.  12  are  the  elastic  cross  sections  evaluated  with  the 
modified  coupling  constants*  -  compare  these  with  the  cross 
sections  from  unmodified  coupling  constants  (solid  curve  in 
fig.  11).  We  find  that  the  cross  sections  are  much  improved 
by  these  changes,  although  they  are  still  not  satisfactory. 

*These  curves  were  computed  directly  from  the 
relativistic  Born  amplitudes  without  partial  wave  decomposition. 
The  Coulomb  scattering  was  included  by  adding  photon-exchange 
diagrams  with  form  factors  as  given  in  (3.1),  (3.2),  (3*3). 
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We  can  thus  understand  the  improvement  in  the  radiative  cross 
sections  in  terms  of  the  corresponding  improvement  in  the 
predicted  elastic  cross  sections. 

Any  model  used  in  NNy  calculation  must  have  the  correct 

soft-photon  behaviour  given  by  Low's  theorem  (the  principle  of 

model  independence) .  In  view  of  the  Burnett-Kroll  theorem  this 

means  that  the  model  must  first  predict  correct  elastic  cross 

sections  and  their  derivatives  (with  respect  to  s  and  t) .  This 

requirement  is  clearly  not  satisfied  by  the  uncorrected  meson- 

pole  model,  as  we  have  seen.  The  correction  used  in  the  present 

calculation  consists  of  using  on-shell  data  to  modify  the  form 

factors  and  to  normalize  the  radiative  cross  sections  as  outlined 

in  section  3.2.  It  should  be  pointed  out  that  the  choice  of  the 

on-shell  point  (as  given  by  (3.1)  and  (3.2))  at  which  we  obtain 

the  normalization  factor  is  arbitrary,  since  there  is  no  unique 

on-shell  point  for  the  expansion  (1.7).  Fortunately  this  is  not 

usually  important  for  two  reasons:  1.  as  long  as  the  choice  of  the 

on-shell  point  (s,t)  satisfies  continuity  requirements  (such  as, 

2  2 

that  ( s , t )  lie  within  the  rectangle  s  e  [  ^  +  P2)  ,  (p3  +  P^)  ], 

t  e  [(p3  -  p1)2,  (p^  -  Pz)2]>  can  affect  I TPPY 1 2  only  in  0(k°)* 
although  it  affects  the  first  two  terms  in  the  expansion  (1.8) 
individually,  it  makes  no  difference  to  their  sum  to  0(k  1) )  [_  28  J- 
2.  we  confine  most  of  our  calculations  to  the  region  where  is 


*It  is  because  of  this  property  that  the  Low  theorem  can  be 
applied  unambiguously  (i.e., 
in  k)  . 


there  is  on  ambiguity  to  two  orders 
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a  fairly  flat  function  of  Hence,  for  soft  photons  the  choice 

of  our  on-shell  point  will  not  affect  results  throughout  the  energy 
range,  but  when  terms  of  0(k  )  can  be  important  we  must  ensure  that 
the  kinetic  energy  corresponding  to  s  =  (p^  +  p^)  is  greater  than 
about  90  MeV. 

2 

A  remark  is  in  order  with  respect  to  the  form  factors,  F  (q  ) 

a 

(table  C.l  with  parameters  given  by  (3.3)).  These  now  carry 
part  of  the  burden  of  correcting  for  the  use  of  the  Born  approx¬ 
imation  in  treating  the  strong  interaction  and  are  much  more 
strongly  damping.  Thus  they  can  no  longer  be  regarded  as  the 
true  strong-interaction  form  factors.  In  the  case  of  the  vector 

mesons  the  cutoffs  are  smaller  than  the  meson  masses,  hence  the 
2 

F  (q  )  are  undefined  on  the  meson  pole.  Note  also  that  the  forms 
v 

2 

of  F^(q  )  used  by  EHM  are  probably  not  realistic  to  begin  with 

even  with  their  parameters,  since  the  values  of  the  ’coupling 

2 

constants’,  g  F  (m  ) ,  are  unrealistically  large, 
v  v  v 

Photon-Emission  Vertex 

With  one  nucleon  leg  off  mass  shell  the  most  general  form  of 
the  vertex  T  ,  diagrammed  in  fig. 14,  contains  two  form  factors 
and  can  be  written  in  the  form  /  44  / 

UCp)  (/=  p  +  k)  -  zzcp)irM  -i 

*  [  f;  cp")  4^  +  ]  o  .i3) 


In  the  limit  as  i'  -*  i  =  m  the  vertex,  T  ,  becomes 


Tm.  ^ 


-2 


(3.14) 
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got 

where  is  the  anomalous  proton  magnetic  moment.  The  form 

factors  in  (3«13)*  therefore,  have  the  mass-shell  limits 


F2.  (  vn'O  7c 
-  'X 


(3.15) 


^  is  a  constant  defined  by  (3.15).  Note  that  the  form  of  the 
nucleon  propagator  given  by  (3.9)  follows  from  the  Ward  identity 


ic  C  p')  S  (  ”  (  p'  ~  HC  p)  (  p'  =  pt  k)  (3-76) 

and  eq.  (3-73)  for  the  form  of  ll  ( p)  f1^  C  p')  . 

Nyman  /P4_7  derived  the  form  factors  from  their  dispersion 

relations  5_7  with  the  assumption  of  threshold  dominance*,  and 

used  them  in  his  soft-photon  approximation  to  investigate  the 

effect  on  the  bremsstrahlung  cross  section.  He  found  that 

whereas  the  9r  distributions  are  little  affected  by  variations 

d  <5 

in  the  form  factors,  the  integrals,  ctJ*i  ,  are  quite 

sensitive.  This  is  not  surprising  since  it  is  known  that 

with  the  static  electromagnetic  vertex,  the  anomalous-moment 

terms  give  a  fairly  flat  photon-angular  distribution,  whereas 

the  charge  terms  are  responsible  for  the  quadrupole  humps 

d  cr 

Furthermore,  the  error  in  due  to  an  incorrect  form 

factor  can  be  much  larger  than  the  discrepancy  between  a  cal¬ 
culation  with  the  static  vertex  and  experiment.  Thus,  until 
the  form  factors  are  determined  from  measurements  it  is  not 
deemed  safe  to  use  those  obtained  from  simple  assumptions  (such 


*The  functions  F 2  - 
pion-product ion  threshold 


(w2) 

i . 


develop  imaginary  parts  above 
e.  when  W_>  m  +  mn  . 


as  threshold  dominance)  in  pptf  calculations.  For  this  reason 
we  use  the  static  vertex  (3-10). 

The  Coulomb  Effect 

The  proper  treatment  would  require  the  inclusion  of  a 
large  number  of  diagrams,  to  order  e3  ,  as  shown  in  fig.  15. 
Furthermore,  the  Coulomb  effect  has  been  studied  by  Signell 
and  Marker  /l8_7  and  by  Heller  and  Rich  /17_/  and  they  found  the 
pure  Coulomb  bremsstrahlung  (corresponding  to  the  first  diagram 
of  fig.  15)  to  be  only  a  small  part  of  the  total  Coulomb  effect. 
Nevertheless  we  included  the  one-photon  exchange  simply  because 
we  already  had  the  machinery  to  handle  vector-boson  exchange. 


CHAPTER  4 


A  CONTRIBUTION  FROM  w  RADIATIVE  DECAY 


4.1  Introduction 

In  this  chapter  we  consider  the  contribution  of  the  uj 
radiative-decay  process  (fig.  16)  to  the  pptf  cross  section. 
This  is  an  internal-radiation  process  and  is  of  0(k)  in  the 
soft-photon  limit.  The  vanishing  of  this  contribution  as  k~*0 
can  be  understood  in  terms  of  the  effect  the  characteristic 
length  of  the  radiating  system  has  on  its  efficiency.  In  the 
case  of  the  external  radiation  the  characteristic  length,  d, 
is  the  distance  the  radiating  nucleon  can  move  with  energy 
imbalance 


Thus  d  —'V'A't  ~  'V  k  and  becomes  increasingly  large  as  k-?>0. 


In  internal  radiation,  on  the  other  hand,  the  photon  is  emitted 
from  within  the  region  of  the  strong  interaction,  thus  the 
characteristic  length  is  a  constant:  it  is  the  range  of  the 
nuclear  force  (or  for  e.g.  the  pion  Compton  wavelength).  There¬ 
fore  the  internal  processes  make  increasingly  less  efficient 
radiators  as  k->0  and  become  negligible  in  this  limit. 

Ueda  /~47 7  had  made  a  rough  estimate  of  the  p  radiative- 
decay  contribution,  and  set  its  upper  bound  at  2.%  of  the  Born 
terms.  We  anticipated  the  (J  radiative-decay  process  to  be 
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more  important  for  two  reasons:  1)  the  ujpp  coupling  is  larger 
than  .p  pp  coupling  (SU(3)  symmetry  with  the  Okubo-Zweig-Iizuka 
(OZI)  rule  and  the  assumption  of  pure  F-type  coupling  imply 
Sujpp  =  3gp^  )  i  2)  the  u>  radiative  width  is  larger  than  the  p 
radiative  width  -  880  -  50  keV  against  35  -  10  keV.  By  this 
argument*  the  interference  from  the  to-p  process  can  be  expected 
to  be  about  15  times  that  due  to  p-/l  process.  If  Ueda’s 
estimate  of  the  upper  limit  is  taken  seriously  then  this  could 
be  an  important  contribution.  Furthermore  the  inclusion  of 
the  process  has  been  found  to  be  necessary  in  calculation 

of  neutral  pion  photoproduction:  it  is  needed  to  reproduce  the 
T  =  -g-  magnetic  dipoles  /P+9_/  • 

In  view  of  the  above  arguments  we  have  carried  out  a  fully 
relativistic  calculation  of  the  process  shown  in  fig.  16.  This 
calculation,  done  in  momentum  space,  is  presented  in  section  4.2. 

A  non-relativist ic  reduction  is  done  in  section  4.3. 

To  ascertain  the  importance  of  the  relativistic  corrections 

at  T,  ,  =  158  MeV  we  have  calculated  the  differential  cross 

lab 

section  for  graph  (a)  of  fig.  16  both  relativist ically  and  in  the 
non-re lati vis tic  limit.  As  the  T  matrix  for  each  graph  is 
individually  gauge  invariant  this  is  a  meaningful  comparison. 

In  section  4.3  we  have  also  derived  the  space  structure 
of  the  non-relat i vist ic  amplitudes,  which  shows  that  the  7T 
propagator  makes  this  a  long-range  process. 


* An  important  factor  excluded  here  is  the  fact  that  the 
caj  has  almost  no  tensor  coupling  to  the  nucleon,  while  the  -p 

has  (  ^  4  )  • 
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4.2  The  Matrix  Element 

We  neglect  the  tensor  coupling  of  the  u;,  which  is  generally- 
assumed  to  be  small  /~50_7*  The  interaction  Lagrangian  is  taken 
to  be 


-  -i  g-n  %r57-TT  \  ~tS^%  CO'1' 


(4.1) 


The  6on T  vertex  is  parametrized  through  the  cJ-?  tt0*-  T  matrix 


K 


UJ7T  ^ 


where  <6^-^  {  )  is  the  tJ  -meson  (photon)  polarization  vector 

and  k  (  pO  is  the  photon  (pion)  momentum.  From  (4.2)  we  get 
the  decay  width  of  the  to  as 


^oo-n-r  C/xwo-p  cu,')  ^  Cy')  ^  P 


(4.2) 


P(lo-^  7i°  r')  ~ 


2 

3 

^  6  7T 


(4.3) 


hV 


UJ 


The  matrix  element  for  the  process  shown  in  fig.  "l6a  is 
T  ~  4  <6  f  fi-')  ^u-3'^yU  44  r)  6:  w  <j  6  p  ^  R  i  ~  Pa'') 


(4.4) 


—  ^  ^  UZ 


(  p3  -  P-Y  -  WVpJ  (  '  pzY-  Yn 

(for  clarity  of  notation  we  do  not  write  the  product  of  the 

coupling  constants  CjTl  <gw  2to7rr  Since 

kp  -  (p ,  -  psY  -4p4-p,Y  (4-5) 

replacing  £  ^  by  Y  causes  (4.4)  to  vanish;  hence  T°^  ±s 


gauge  invariant . 

Introducing  the  following  notation: 


45 


^  (  L  J  "  [  (  P  L  ~  Pj')  ^2uol 

■^Ti  (lj -  r  c pt  -  pj)  -  h^\j 


^V6  p 


we  can  write 


\ 

>  (4.6) 

) 


I  T)coCi  3)  bjr  ^5  ki'l  (4.7) 


Similarly 

^  o')  -fap  ^l.^)  I)ui  ("2,  SO  ZD7T  0,  Tj-  66,^)  (4.8) 


where  ^  6p^p  feVp3-p^  . 

The  contributions  of  the  exchange  graphs,  fig.  16c  and  l6d,  are 
obtained  by  the  interchange  of  labels  3^4  in  the  expressions 
(4.7)  and  (4.8).  The  complete  matrix  element  is  then 


tl 


— —  c  1  <0) 


+• 


-r  ClO 


-  T 


c»0 


T  CM') 


(4.9) 


This  was  evaluated  numerically  and  added  to  the  OBE  matrix 
element  as  given  in  section  3*2.  The  rest  of  the  calculation 
is  the  same  as  in  the  OBE  case. 

zz  1  ^  TiS  r, 

-  although  tedious  -  was  also  done  analytically.  The  two 
methods  gave  identical  results. 

The  coupling  constants  used  in  the  evaluation  of  T^n^.  were 


As  a  check,  the  calculation  of  ■=f  4. 

Spins  6 


9* 

H  71 


^  to  —  5.4  ^lotty  -  H.<?7  t  /  0  2  (4.10) 

TU  71 


IH.  O 
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4.3  Non-Relativist ic  Reduction  and  the  Coordinate- 

Space  Structure  of  the  Matrix  Element 

The  non-relativist ic  reduction  of  the  matrix  elements 
(4.7)  and  (4.8)  is  done  by  means  of  the  approximations 


ui  <r^Uj — *  UL  y„  Uj  * 

[  Cpi  ~  pi)  - y»' 4]  *  -["£2  +  »4]  =  <p»- p,Y 


(4.11) 


[  ( p4- *  [■p%^17r]  pM  =  (pl1-p^y 


The  above  replacements  yield 


T 


(  I 


Y"> ,  v. 


“  I 


& 

Q  yy\ 


(kff) 


( 1  +•  ^4)  (~p 2  +  rr  ") 


(4.12) 


r 


ciO 


n.  Y\ 


4 


2  m 


■0? 


*?) 


+■  h^lo')  +  v4u>^) 


(4.13) 


Each  of  these  terms  represents  a  conserved  current.  In  the 

(la) 

arguments  that  follow  we  consider  only  the  term  T  .  To 
obtain  the  non-relativist ic  cross  section  we  use  the  identities 


and 


L  =  /I*  -  -T)x 

L  \(  xf  'Z-a  Yi)!1  = 

SpinS 


(4.14) 


for  any  3-vector  A. 


We  get 


(4.15) 


4  7 


■L  1 1  C 


Spiv>5 


2  2  \2. 
+  h^TT  ) 


The  relativistic  version  of  (4.15),  derived  from  (4.7),  is 

Z  \  t0,Y  =  J-u  [( k- ?)((?*■  k)(P,-p) 


H 


SpinS 

£■* 


2  w 


-h 


(Pi  ^)(p3-  p))  ~  P2(pr  kXp2- k)  ~  (k-p)*] 

*  L  Cpw*K)-^ J  f  Dw0.3)  DjtCs^l 


(4.16) 


Eq.  (4.15)  can  also  be  derived  by  non-relativist ically  reducing 

(4. 16)  .  The  cross  section  is  obtained  by  multiplying  (4.15)  or 

2 

(4.16)  by  (  ^ 7i  3to7tr  )  x  (phase-space  factor). 

In  table  4.1  we  show  both  the  relativistic  and  non- 

relativistic  values  of  the  differential  cross  section  as  obtained 
(la) 

from  T  above.  It  is  clear  that  at  T.  ,  =  158  MeV  the 

lab 

relativistic  corrections  are  important 


©  (deg) 

Relativistic 

Non-relat ivist ic 

0 

0.54 

1.14 

30 

1.89 

3.02 

60 

1.83 

2.88 

90 

0.83 

1.52 

120 

0.16 

O.38 

150 

0.27 

0.00 

180 

1.03 

0.29 

Table  4.1.  dcr/djl^  d  (9r  (jjl  b/s  r1-  ra.d')  <3^6  "i5ab  /£><?  A/eV 

(e^eH--i  0°) 
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Coordinate-Space  Structure 

Define  the  coordinates 


~r,  ,  R  -  i.  ^  -*-"*0 


Q  =•  i  (^p,  -pO  ,  p  “  ~p>  +-  f 


(4.17) 


and  analogous  expressions  for  the  primed  quantities  which  refer 
to  the  final  state  coordinates  and  momenta.  The  Fourier  trans¬ 
form  is  then  defined  as  /51_7 

T(r,  R,"?'  R')  =  fay,  $  -€ 

*  §C3V  P  -  P'-t)  € 


i P • R  -7 


oTP  d'Q  cJl4P'aV 


iP-R-iQ-?  l3r,  fi  a/  |3  .  (4.18) 


(la)  -** 

Inserting  T  for  T(P,Q;P,,QI)  from  (4.12)  one  gets 

n  •  r  • 


—r-  ( 1  tt~) 


0?,  t  J r‘,  R')  =  ^  ^  R  S C^(R-R')  (4.19) 

*  $c,Y?-r)  tc'*V) 


where 


T^CA  =  -£  •  ClTvf^  ♦  *;?)] 


where 


-e 


-a-'? 


1  =  )  [(T*-ktY+«,;][(7-nrT^t] 


(4.20) 


(4.21) 


For  low-energy  photons  (k-*0) 

=  -£l_-  f^ax^.vjl^ 

£2  im 


where 


x  _ _  /  JL  -  -£ 

XCrjO )  -  Hn  (  ^X-^7r) 


v~ 


v~ 


(4.22) 


(4.23) 


/1bx 

The  coordinate  space  expression  for  T'1  ;  is  obtained  by  the  simpl 

replacement  CT,  — =*  CT^  .  Clearly  the  amplitude  has  a  short-range 
part  and  a  long-range  tail.  At  r  =  0,  I(r,  0)  =  1/(4tt)  (m^  +  m^  ). 
Thus  the  large  mass  of  UJ  does  indeed  suppress  the  overall  con¬ 
tribution  of  this  internal-radiation  process. 

4.4  A  Discussion  on  the  Coupling  Constants 

First,  we  point  out  that  in  the  calculation  of  the  OBE  con¬ 
tribution  we  were  constrained  to  using  the  parameters  of  the  EHM 
model  (albeit  with  some  modifications  as  discussed  previously)  in 
order  to  not  disturb  the  fit  to  the  elastic  data.  This  was  not 
the  case  in  the  other  contributions  we  considered:  here  we  were 
free  to  use  the  known  neutral-meson  masses,  and  coupling  constants 
consistent  with  physical  principles  and  experiments.  For  example 

we  disagree  with  the  constraint  imposed  by  EHM  on  their  vector 

2  2  2  9 

meson-nucleon  coupling  constants  that  g^  :  g^  :  g^  =  1:  —  :  9. 

We  do  not  expect  the  4>  ,  which  is  a  pure  13  state  in  the  ideal 
nonet-mixing  scheme,  to  couple  to  protons,  which  have  no  strange 
quark  content.  This  is  the  statement  of  the  OZI  rule  /  52_7* 

We  shall  now  examine  the  constraints  imposed  on  the  coupling  con¬ 
stants  by  vector-meson  dominance  (VMD)  and  unitary  symmetry. 

Consider  the  coupling  of  the  nucleon  to  the  electromagnetic 
field  as  in  fig.  17a.  The  invariant  amplitude  (from  the  Dirac 
coupling  alone)  for  this  diagram  can  be  generally  written  as 


NCpO  ~(M  «  [  f r,s  (0)  -  t3  *7  (o')]  NCp) 


(4.24) 


where  F^(q2)  and  F^(q2)  are  the  isoscalar  and  isovector  charge 
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form  factors  (F^(0)  =  F^(0)  =  -5-) .  If  we  assume  the  isoscalar 

part  of  fig.  17a  to  be  dominated  by  to  and  exchange,  and  the 
isovector  by  p  exchange  (fig.  17b)  the  amplitude  is 


6?  V 


M. 


<£  wn 


u> 


-fii? _  +  -3^  +r  f 

P  P  V  ^Tp  ^p-^J^rO 


(4.25) 

NCp) 


■€  ^ , 


where  ^  Tv  (v  -  60_,  <j}  p")  is  the  V-T  coupling  constant.  Com¬ 

parison  of  (4.24-)  and  (4.25)  gives  the  constraint 

3^ 


H- 


V, 


U> 


3^  —  1 

Yf 


)  (4.26) 


3L 


1 


thus  VMD  makes  no  statement  on  g w  or  g^  separately.  We  shall 
now  look  at  the  SU(3)  constraints. 

The  SU(3)  structure  of  the  electromagnetic  current  in  fig. 


17  is 


I«.m.  -  Y  b  L  ( $  3  £ 


r~-  I  n,j  y  -'J  (4.27) 

‘•jJ  V  5 

(  are  the  SU(3)  antisymmetric  structure  constants)  since  it 

must  transform  like  Q  =  F3  +  Jf  Fff  .  Comparison  of  (4.27) 

with  that  due  to  the  general  SU(3)  invariant  F-type  BBV  coupling, 


£  8^  e,  v, 

Cjj 

to  the  combination 


,  shows  that  the  photon  couples  only 


P 


JT 


(4.28) 


where  CJy(u£)is  the  octet  (singlet)  isoscalar, 
ideal  mixing 


If  we  assume 
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\ 

>  (4.29) 


the  combination  (4.28)  is 


P' 


jT  (f> 

3 


+  3  " 


(4.30) 


If  we  ignore  the  mass  splittings  (recall  the  Vy  coupling  constant 
is  €  / (2  '&VS)  )  the  three  independent 

duced  to  one  by  the  constraint 


’  oo 


are  re- 


V  *  X 


2. 


=  "T  :  2.  :  l 


(4.31) 


Eqs •  (4.31)  and  (4.26)  can  now  be  combined  to  give 


^  J~2~  “  3  Jp  (4.32) 

The  OZI  rule  (^^  =0)  then  gives 

4e  :  3l,  :  3-  =  '  :  0  '•  1  (4'35) 

This  constraint  can  be  generalized  to  allow  for  both  f- 
and  d-  type  coupling  of  the  meson  octet  to  the  baryon  octet  as 
follows.  The  general  SU(3)  invariant  coupling  of 

the  meson  nonet  to  the  baryons  is 

J°  -  L  °\f  t  ^  B i(i  v'jBj  +('/-•«') 

*  _  (4.34) 

19.  B. 


x  Si  (i-  5/.1  Bj  J  +  L  g0  B,  isv)  8£  \/ 

i-i 


o 


^  ^-j=~ )  P  t5  ^  3°  PP^o  +■  other  terms 
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With  the  mixing  as  given  by  (4.32)  this  is 


£  =*  C  g  P  J5  JIJ  PPrf 

■*-  C  gp  +•  9°  f|  1  P  Pu 

3  -* 

We  can  read  off  g^  and  g^  from  (4.35);  setting  g^ 


(4.35) 


0  gives 


9®  "  9p  /f  ("3  -^<0 


(4.36) 


which  reduces  to  (4.33)  when  c<  =  0.  Note  that  in  this  argument 
we  did  not  use  VMD. 

2 

Our  choice  of  g‘  LU  (4.10)  is  dictated  by  the  reasoning  pre¬ 
sented  above;  fj^oTiy  was  determined  from  (4.3)  with  r(uj-^7nf) 
=  0.89  MeV  /"53  7. 


CHAPTER  5 


RESULTS  AND  DISCUSSION 


The  coplanar,  symmetric  cross  sections,  ^ ^/d Jli  dSlc,  d&r } 


m 


at  various  lab  energies  and  angles  0^,  0^  are  shown  plotted  i 
figs.  18-34.  The  curves  marked  "OBE"  are  the  results  of  the 
OBE-external-emission  calculation  outlined  in  chapter  3«  We 
were  not  able  to  fix  the  sign  of  the  cj-7T°  amplitude  (eq.  4.9) 
relative  to  the  OBE,  we  therefore  show  the  results  for  both 
signs.  These  are  marked  "OBE  -  internal”  in  the  figures.  The 
cross  sections  were  calculated  at  10°  intervals  in  9*  ;  the 
curves  were  then  drawn  by  interpolation*. 

Fig.  19  shows  a  comparison  with  the  preliminary  results  of 
the  Triumf  experiment  of  Beveridge  et  al  /10_/  and  two  other 
calculations:  a  Hamada-Johnston  potential  calculation  of 

Bohannon  (marked  ”HJ”)  and  a  soft-photon  approximation  calculation 
of  Fearing  (marked  "SPA").  The  experimental  and  the  two  sets  of 
calculated  points  were  adapted  from  fig.  6  of  ref.  31**.  The 
differences  among  the  calculations  are  still  too  small  relative 
to  the  experimental  errors.  The  SPA  is  almost  consistently 
low,  while  the  HJ  calculation  is  high;  however  part  of  this 

£  £  0  j*  0  q  probably  due  to  relativistic  corrections  not  model 


*Since  all  observables  are  generally  differentiable  fun¬ 
ctions  in  0y  they  must  have  a  zero  slope  at  0*  =  0  and  l8o  in 
the  symmetric  geometry.  As  can  be  seen  in  the  figures  this 
requirement  was  sometimes  not  met  by  the  interpolation  procedure. 

**The  OBE  calculation  shown  in  this  figure  did  not  include 
the  correction  described  in  chapter  3» 
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dependence.  Relativistic  corrections,  as  discussed  by  Liou  and 


Sobel  /16J7»  lower  the  cross  section  and  are  largest  at  forward 


and  backward  angles.  A  comparison  of  the  OBE  and  the  SPA  cal¬ 
culation  provides  a  better  indicator  of  model  dependence  since 
both  are  fully  relativistic.  It  appears  unlikely  a  definite 
conclusion  regarding  potentials  will  be  made  on  the  basis  of 
this  experiment.  Fig.  37  is  a  plot  of  the  off-shell  parameters 


and  the  photon  energy,  k,  as  a  function  of  ©v  .  It  shows  that 
the  virtual  nucleon  is  farther  off-mass  shell  at  the  forward  and 
backward  angles  -  the  region  where  the  largest  model  dependence 


is  observed.  On  the  other  hand  the  photon-energy  distribution, 


k(©v  ),  does  not  appear  to  be  a  useful  indicator  of  where  (in 
©  *  )  model  dependence  should  be  observed.  The  problem  here  is 
that  we  do  not  know  how  the  higher  order  terms  in  the  k  expansion 
depend  on  ©^  to  be  able  to  predict  their  importance  -  even  when 
we  know  the  photon  energy. 

The  co  radiative-decay  contribution  (labelled  "internal”  in 
fig.  19)  is  not  negligible  for  this  experiment:  it  is  largest 
in  the  forward  and  backward  photon  directions  (about  10%  at  the 
forward  direction).  We  are  thus  faced  with  having  to  include 
meson-current  correction  as  well  as  rescattering,  Coulomb  effect 
and  relativistic  corrections  before  meaningful  comparisons  can  be 
made  between  potential-model  calculations  and  the  experiment. 

Should  experiments  be  performed  at  larger  energies  and 
smaller  proton  angles  the  meson— current  corrections  will  be  even 
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larger.  At  T1&b  =  400  MeV,  0^  =  =  12°  (fig.  15)  the  ^-rr° 

contribution  is  about  27%  in  the  forward  photon  direction  and  25% 
in  the  backward  photon  direction;  the  effect  on  the  integrated 
cross  section  d«/dSllCLJZH  is  18%,  as  can  be  seen  in  table  5 • 1 • 
Fig.  29  shows  that  at  158  MeV  ©^  =  ©^  =  30°  the  lo-7t°  con¬ 
tribution  is  a  fraction  of  the  size  of  the  error  bars;  for 

©2.35°  =  158  MeV)  it  is  negligible.  This  behaviour  is 

expected  since  the  photon  energy,  k,  decreases  with  increasing 
proton  angle  ©  and  the  amplitude  is  of  0(k)  in  the  soft- 

photon  limit. 

The  99  MeV  calculations  (figs.  3‘1  and  33  and  table  5-1) 

are  actually  beyond  the  range  of  this  model  and  should  not  be 

taken  seriously.  The  maximum  values  of  T^ab,  as  computed  from 

s  =  -J  /Tp1  +  P2)2  +  +  P4)f7,  are  71*3»  63.6  and  58.4  MeV  for 

©  =  35°,  30°  and  25°  respectively.  Recall  that  the  present  OBE 

model  is  not  valid  when  T..  ,  < 100  MeV.  What  the  calculations 

lab 

of  figs.  14  and  16  show  is  that  the  one-photon-exchange  con¬ 
tribution  is  indeed  negligible  as  observed  previously  by  Heller 
and  Rich  /17_7  and  Signell  and  Marker  /1<3_/. 

In  figs.  21  to  28  we  present  a  study  of  the  various  con¬ 
tributions  to  the  200  MeV  ©  =  1 6 .3°  calculation.  Figs.  21  and 
22  show  how  each  contribution  affects  the  coherent  sum,  while 
fig.  23-28  show  the  distributions  due  to  each  contribution 
individually.  We  observe  that  the  scalar-meson  exchange 
enhances  the  cross  section,  particularly  at  the  peaks  of  the 
quadrupole  humps,  while  the  vector  mesons  reduce  the  cross 
section.  The  general  enhancement  at  the  backward  photon 
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direction  is  due  to  the  phase-space  factor,  which  has  a  minimum 
at  9r^  70°. 

We  get  a  better  understanding  of  the  shape  of  these  dis¬ 
tributions  from  fig,  34,  which  shows  the  contributions  from  the 
Dirac  (M^-p  =  0M  curve)  and  anomalous-magnetic  moment  (M  ocP  term” 
curve)  coupling  of  the  photon  to  the  radiating  proton.  The 
Dirac  term  contains  a  charge  and  magnetic-moment  part  while  the 
anomalous  term  contains  only  a  magnet ic -moment  part.  We  find 
that  the  main  contribution  is  due  to  the  magnetic-moment 
amplitude,  which  is  roughly  isotropic.  The  charge  amplitude 
still  gives  rise  to  the  quadrupole  distribution,  but  this  is  not 
large  enough  to  be  seen  in  the  total  distribution. 

Finally,  we  come  to  the  comparison  of  our  calculation  with 
the  results  from  Orsay  ^34/.  This  is  one  of  the  two  experiments 
where  cross-section  data  were  obtained  at  small  and  unequal  proton 
exit  angles  and  9^.  The  discrepancy  with  previous  calculations 


has  been  discussed  by  a  number  of  authors  /8 ,  17,  55,  56/.  It 
has  been  suggested  (by  Jovanovich  /17_/)  that  this  experiment 
possibly  indicates  a  failure  of  the  Hamada-Johnston  potential 
off— energy  shell.  However  our  calculation  (figs.  35  and  36) 
was  done  in  a  completely  different  way  and  yet  we  find  essentially 
the  same  discrepancy.  Furthermore,  we  find  that  this  cannot  be 
accounted  for  by  the  u^-m0  current  correction  or  relativistic 
corrections*.  Jovanovich  /  1 7__/  has  suggested  that  the  42  MeV 


* Relat i vist ic  corrections  in  the  cont_ext_of  this  experiment 
have  also  been  discussed  by  Celenza  et  al  /55_/* 
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Manitoba  experiment  shows  the  same  type  of  discrepancy  at  small 
and  unequal  proton  angles.  The  Orsay  results  thus  remain  as  an 
enigma. 


Present  work 


Hamada- 


T 

lab 

(MeV) 

93’  9 4 
(deg) 

OBE 

OBE 

+intemal 

OBE 

-internal 

SPA 

[27] 

Johnston 

Potential 

[191 

Experiment 

99 

25 

4.19 

3.77 

±  0.23**  [48b] 

30 

5.94 

25.0* 

5.14 

±  0.22**  [48b] 

35 

8.96 

*7.3* 

9.0 

±  0.3**  [48b] 

40 

19.8 

*17* 

18.8 

±  1.2**  [48b] 

156 

15,  18 

4.05 

3.96 

4.17 

15,  21 

4.28 

4.21 

4.40 

15,  24 

4.40 

4.33 

4.50 

15,  27 

4.41 

4.35 

4.50 

158 

30 

8.71 

8.65 

8.84 

7.9 

9.15 

10.2 

±  1.7  [4c] 

7.9 

+  1.6  [4b] 

35 

12.7 

11.8 

13.2 

14.7 

+  2.5  [4c] 

12.4 

±  2.5  [4b] 

40 

32.5 

33.0 

23.8 

±4.8  [4b] 

200 

12 

3.42 

3.25 

3.66 

16.3 

4.85 

4.66 

5.13 

25 

8.38 

8.21  . 

8.69 

30 

11.6 

10.0 

12.1 

13.0 

±  2.4  [6] 

35 

17.1 

15.3 

17.5 

14.0 

±  2.7  [6] 

40 

46.4 

46.3 

46.5 

49.1 

29.0 

+6.0  [6] 

400 

12 

10.7 

9.09 

13.1 

Table  5.1  Integrated  cross  sections  da  /dSL%  d  (yub/srl>). 

"OBE"  is  the  present  one-boson-exchange  external-emission  calculation. 

"OBE  ±  internal"  includes  the  radiative-decay  contribution 

(chapter  4),  coherently  added  with  either  relative  sign.  "SPA"  is  the  soft- 

photon-approximation  calculation  of  Nyman  [27] . 

*  Estimated  from  figs.  13  and  14  of  ref.  19. 

**  Additional  normalization  factor  of  1  +  0.02  to  be  applied. 


CHAPTER  6 


A  CONTRIBUTION  FROM  A  EXCITATION 


6.1  Introduction 

In  this  chapter  we  consider  the  simplest  radiative  pro¬ 
cesses  involving  A  excitation  with  tt  and  p  exchange  (fig.  38). 
This  contribution  is  of  particular  interest  at  the  larger  incident 
energy  of  730  MeV  of  the  UCLA  experiment  where  in  certain 

kinematic  situations  it  can  dominate  the  cross  section.  Recently, 
Tiator,  Weber,  and  Drechsel  /_€>3 J  published  the  results  of  a 
calculation  in  which  the  above  mentioned  A  graphs  were  added 
incoherently  to  the  first  term  of  the  Low  expansion.  In  the 
present  work  these  A  graphs  are  added  coherently  to  the  OBE 
model  of  chapter  3  and  the  results  are  compared  with  the  UCLA 
data.  We  find  the  interference  terms  to  be  very  important,  and 
when  they  are  included  some  of  the  conclusions  made  in  /6j^/ 
need  to  be  modified. 

To  estimate,  a  priori,  the  importance  of  this  process  we 

evaluated  the  f our-momenta-squared  of  the  virtual  baryon: 

(Pl  -  k)2,  (p2  -  k)2,  (Pj  +  k)2,  (Pif  +  k)2.  At  Tlab  =  200  MeV, 

Q  =  ©,  =  12°,  the  smallest  and  largest  value  of  these  invariants 

3  4  _ 

is  0.64  GeV2  and  1.06  GeV  .  The  on-shell  value  is  m  =  0.88 
GeV2.  Compared  with  m2  =  1-52  GeV2  the  off-shell  excursions 
are  not  large,  therefore  we  expect  the  nucleon-pole  terms  to 
dominate.  On  the  other  hand  at  Tlab  =  730  MeV;  9^  =  67  , 
ip  =  -1.09°,  6*  =  50. 5°,  f- r  =  180°(G7  geometry  of  the  UCLA 
experiment  /l1_7)  the  smallest  and  largest  value  of  the 
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2  2 

invariants  is  .52  GeV  and  1.27  GeV  .  The  largest  value  is  now 

2 

somewhat  closer  to  mA.  We  therefore  expected  that  the  A  con¬ 
tribution  may  be  significant  at  the  energy  of  the  UCLA  experiment. 

6.2  The  Matrix  Element 
ANn  Vertex  and  A  Propagator 

The  most  general  charge  independent  ANt\  interaction  can 

be  written  /*57_7 


=  g *  t*  T  &  ±k.c. 

- 

where  ^  +  f  k  (l  +  ^  Z)  /4  +  ZJ  T^X  , 

A,  Z  are  complex  constants, 


(6.1) 


T  are  4x2  matrices  in  isospin  space  analogous  to 
the  't  matrices  in  NNn  coupling, 


"  .+■+•  1 

b 

II 

A 

Av 

and  aL  _ 

'N 

? 

A° 

n 

A" 

The  explicit  isospin  structure  can  be  obtained  by  forming  the 


invariant 


The 

TTm 


L  \  a*)  <  /  4 

t*  t 

<f>  (m) ,  m  =  +1,  0,  -1 

;  in  terms  of  the  c 


i  £*>  %  (t)  HA 


,  are  annihilation  operators  of 
artesian  fields  ,  they  are 


(6.2) 


the 


~  ACm'')  •  <j> 

£(+\)  ~  Jz'  (~\j  L)  ^ 

6(A~- 

6  (o)  -  (0,0,  l) 


The  matrices  can  now  be  read  off  from  (6.2)  (up  to  an  overall 
real  constant*): 
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I  L  vyi  -L  I  4^  (vyC) 


*\  t 


£ 

h 

]/2 

->4 

~34 


£60 
O 


JT 


0 

^7? 

~£r(6)  Jf 

£(-l) 


(6.4) 


The  coupling  (6.1)  is  commonly  and  conveniently  written  in  terms 
of  doublets  N*i  =  %  Ti  (or  N*.  =  ).  According  to  (6.4) 
they  are 

,v* 


M 


A+"  -  AV/T 


\A/f3  -  A' 


K) 


M 


(6.5) 


3  “ 


~li 


at 

A° 


The  general  spin  3/2  propagator  is  a  sum  of  the  Rarita- 
Schwinger  propagator  /  72_/  and  a  non-pole  term  that  depends  on 
the  parameter  A,  defined  above  /  5 7_/ ;  this  term  vanishes  when 
A  =  -1 .  The  total  Lagrangian  (free-field  part  plus  interaction) 
is  found  to  be  invariant  under  the  point  transformation 


♦The  minus  sign  in  (6.4)  with  the  Condon  and  Shortly  sign 
convention  used  here  results  in  the  doublets,  Nf  (6.5),  as  given 
in  the  literature. 
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A  =  +  «-  VMh  ^4 

a'  =  f/i  -  + 

a  =  ^  (6 

<fj,  =  ^7T 

The  S  matrix  is  therefore  independent  of  A;  without  loss  of 
generality  we  took  A  =  -1  to  get  the  simplest  A  propagator: 


5pw  (fp*') 


_  j_  /*  +  -^pw  - ^  A-jSjo  p’ v'APp 


P  “  >™A 


(6.7) 


+■  — a  P*P  pv  ^ 


To  take  the  width  of  the  A  into  account  and  to  remove  the 
2  2 

p*  =  mA  singularity  we  made  the  replacement  in  the  above 
propagator: 


XI  X 

P  - 


-  L  r  #0*0 


(6.8) 


where 


-P  ( pO  =  I?" 


*1 


A 


m 


\  -  A 


p 

/\  =  (threshold  value  of  p*  for  the  decay 

2 

A~ *  Nn)  =  (m  +  m^  ) 

I  '  is  the  A  width  =  115  MeV 


Olsson,  Turner  and  Osypowski  /  57/  found  from  their 
analysis  of  tv  N  scattering  that  Z  must  lie  in  the  range 
-0.8-  Z^O.  Olsson  and  Osypowski  /~58_7  further  restricted  Z 
by  comparing  with  a  multipole  analysis  of  pion  photoproduction 
in  the  A -resonance  region:  they  concluded  that  Z  =  0  -  ■£. 

We  have  nevertheless  used  Z  =  — J  since,  along  with  A  =  -1 ,  it 
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considerably  simplifies  the  AN  7T  interaction;  i.e.  we  have 
used 


=  9*  (  F*  M  \TTi  +  w  X  TTi 


For  the  case  of  A  p  Tt  coupling  this  reads 


(6.9) 


(6.10) 


where  <^>rt  are  the  A+  ,  proton,  and  TT°  fields 

respectively. 

We  can  now  write  the  matrix  elements  corresponding  to 


diagrams  I  and  II  of  fig.  39: 


^71  C^)  = 


F7~ 


A 


\f £ 

thus  ~dM  (j)T[  LoF)—>  -l  ^ 


+-  b* 


in  both  diagrams  and 


lr  ~  “  3*  i|  ^P^ 

Ttt  -  -3"  if  -(P^  C-c^  — fpO 


(6.11 ) 


AN  *  Vertex 

The  general,  gauge-invariant  ANY  couplings  were  investigated 
by  Gourdin  and  Salin  /"58_7.  They  showed  that  for  free  fields 
there  are  two  such  couplings: 


and 


f-  +  ^  f 


V  1  w 

the  first  of  these  producing  mainly  magnetic  transitions,  the 
other  primarily  electric  transitions.  From  analysis  of  pion 
photoproduction  experiments  it  is  known  that  the  process 


64 


^  P  ^  N  tt  is  dominated  by  the  magnetic  transition*. 

We  consider  only  the  dominant  coupling.  Olsson  and 

Osypowski  /~£>0_/  add  an  off-shell  term  to  this  coupling  so  that 
the  total  Lagrangian  is  invariant  under  the  point  transformation 
defined  in  (6.6).  (This  extra  term  vanishes  when  the  A  field 
is  on  mass  shell).  The  resultant  interaction  is 


of 


Ay  m 


-  -  e.  C  AJ 


#3 


M- 


o 


r- 


N  F 


In.  C. 


(6.12) 


where 


[Y  i(l  +  vr)fl 


jY^Y, 


Y  is  an  arbitrary  constant ; 

and  A  is  the  constant  defined  before,  which  we  set  as  A  =  -1 . 

The  authors  determined,  from  a  comparison  with  pion-phot oproduc t ion 
multipole  analysis,  that  Y  =  %  -  t|;.  For  simplicity  we  took 
Y  =  --g-,  the  value  that  along  with  A  =  -1  removes  the  off-shell 
term  in  (6.12).  For  the  p  (6.12)  reads 


A 


Up 


-  € 


C 


(6.13) 


—  —  H* 

(In  the  notation  of  Peccei  /6l_/  C  =  ,  where  m  =  nucleon  mass). 

We  can  now  obtain  the  matrix  elements  for  diagrams  I  and  II  of 


*For  example,  the  analysis  of  Berends  and  Donnachie  /_59_J 
suggests  that  the  isospin  3/2,  electric-quadrupole  amplitude 
vanishes  near  the  resonance  region.  Note  that  because  of  parity 
and  angular-momentum  conservation  only  two  multipoles  can  con¬ 
tribute  to  r +  A— >  tt-i-a)  .  magnetic  dipole  ,  and 

electric  quadrupole  (the  superscript  is  twice  the  total 

isospin) . 
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fig.  40: 


^  (  a  4  a+ 


ftp.  ~  L  ku  6  (-a.  -t 


-f  le  ■ 


-t-  c 4r  -e 


i.  k: 


for  both  diagrams  ?)p/) 


M 


^  L  6m 


L  ~  ~*CJt  -V^Lts  u^Cp*')  i(kp4w-kv^ 

-  £  C  Jj  ZL  CpZ  ( <£  -  it  </)  ys  U-M(p*) 

(6.14) 

7~ir  ~  ~  C  £$  u'*(p*')  (  £  ~X  Tr  a  fp} 


A  Sip  Vertex 

In  analogy  to  the  Ashf  coupling  (eq.  ("6.12)  without  the  off- 
shell  term)  we  took  the  A  Si  p  coupling  to  be 

Amp  -  ~C'  +  4.  «r.  (6.15) 

In  particular 

A+Pf*  =  /!  c'[ZM^trs^  -^y-irs^] 

*  fa?*  -  (6-l6) 

We  can  relate  the  coupling  constant  C'  to  C  by  assuming  vector- 
meson  dominance  (VMD). 

* 

Consider  the  coupling  of  a  photon  to  the  A* p>  under  the 
assumption  of  VMD  (fig.  41 ) ;  note  that  isospin  conservation 
forbids  the  to  and  ^  ,  so  that  only  the  p°  is  involved, 
matrix  element  for  this  process,  according  to  (6.16),  is 


The 


66 


T  = 


-Z 

•e  m  p 

V 


Cv'i  ?\ 


*  u^fpO  C<£ff^  -  K  £(f-)^)rs  uCp)  (6.1?) 


making  the  replacement 


T.  e 


3  •> 


fecp> 


with  k  =  0  we  have 


Cp^A  ^(p'W 


3* 


I?  1? 


Y* 


(6.18) 


7"  =  -Jl£L  JT  UCM(p^  (  kyU  ~  ^  6yS)  YS  U-Cfi) 

0.  Vp  V3 


(6.19) 


Comparing  (6.19)  with  7Jj  of  (6.14),  and  using  Yp  =  cjp  (eq. 
4.26) we  get 

C/  =  2gf,<T 


(6.20) 


The  Bremsstrahlung  Matrix  Element 

Putting  the  various  pieces  together  we  get  the  following 
matrix  element  for  the  7T -exchange  processes  shown  in  fig.  42. 


'OO  + 


=  -2 


+■  c/O 


^  €  C  £)jt  F(  ~  Spw^pa^^  % 


-  Z  ,  1 

$  ~ 


(p3+i)  -  vn\  +-  L  YY)&  p  ffol+h')  ) 


(6.21) 


-  S”pu  (fp,-k)(V*  3^  Ys_  u,  Ys  u^) 

(pr^)Z  $  C(pr^f) 

where 

SpvCp9  =  C^+m^f-^gpw  -t  -UCTrppt,-Tr„p*p')  +  2  p'p  pt 


is  defined  in  (6.8) 
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2 

F(q  )  is  a  form  factor  appearing  at  each  vertex  involving  an 
exchanged  meson.  Following  [  64  J  and  [62,  J  we  take  the 
monopole  form  (the  same  for  each  vertex) 


F<yy  =  az/(az  -  y) 


(6.22) 


with  A  =  1.5  GeV^ . 


For  the  p -exchange  processes  the  corresponding  matrix  element  is 
1  (!*')+  ci  kA 


I  A,. /•.i.'.  =  "^g,uf>Vs 


(  Pj4  fey  -  mi  *  <■  mA  P  ^(1  p3*-fey) 

+  (A  ^  ~  A  jA)  ^  5pv<fP ,  -  feu  - 

Cp.-feV  ~  A  +  t  mA  P  JA[W) 

x  ■%  J  u,  uH  (Vg  - 1  a-^  -Vjp  \  u4 

5  ^  ' 


(6.23) 


<6 


Values  of  the  coupling  constants: 


WTi 


IH 


Ip  -  <?.  6 

Hr, 


-3.7 


z  ~5 

[\  (ANtt)  ~  I.  875  *  JO  A/eV  from  the  A  decay 

1/  7T  width  (appendix  F). 

-  0,35%  myT'1  (with  ^  137  Me\0 

.  ,  x  N  _  U  a  7  //o  \  from  the  multipole  analysis  of 

I  C  |  (/VAT)  -  pion  photoproduction  of  Olsson  [  66  J . 

(WAp)  -  0.  C  <jf>  from  the  VMD  argument  given  above  (6.20). 

An  important  aspect  was  the  determination  of  the  relative 

signs  of  these  constants,  as  the  interference  terms  critically 
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depend  on  these.  The  relative  signs  are  as  follows:  g*gr,  <  0  - 
determined  from  the  quark-model  constraint  (E.  19);  C<0  -  from 

the  quark-model  constraint  (G.5);  gp  >  0  (£w  r - Vi* f 

+  derivative  coupling)  -  from  analogy  with  ;  C'  <  0 

from  the  VMD  constraint  (eq.  6.20).  These  sign  relationships  are 
summarized  in  fig.  43. 

6.3  Results  and  Discussion 


k  are 


The  differential  cross  sections 


plotted  in  figs.  44-49;  the  results  are  compared  with  data  from 
the  five  photon  detectors,  G7,  G8,  G10,  G11,  and  G12  of  the 
730  MeV  UCLA  experiment  /  11_/.  As  expected,  for  small  photon 
energies,  k,  the  predicted  results  show  a  1/k  dependence  and 
are  in  approximate  agreement  with  the  EED  calculation  of  Nefkens 
et  al  /II _7* •  At  larger  photon  energies  our  cross  sections 

increase  with  the  photon  energy  thus  showing  a  minimum  at  some 
intermediate  k.  This  behaviour  is  qualitatively  in  agreement 
with  that  of  the  complete  soft-photon  calculation  of  Fearing  /  3 \7 • 
It  should  be  pointed  out  that  the  presence  of  higher  order  terms 
in  the  expansion  in  k  of  the  cross  section  is  not  necessarily  due 
to  off-shell  effects.  These  may  be  due  to  the  expansion  of  the 
phase-space  factor  or  the  momentum  variables  in  the  amplitudes. 
Indeed,  the  phase-space  factor  is  singular  at  the  maximum 
kinematically  allowed  k. 

The  above  applies  to  the  OBE  calculation  with  and  without 

♦The  ’external-emission  dominance'  or  'EED'  is  a  calculation 
from  the  leading  term  of  Low's  expansion. 


the  A  term  added.  We  turn  our  attention  to  the  region  of  larger 
photon  energy  (k  80  MeV),  where  the  A  terms  become  significant. 
We  can  write  the  total  cross  section  as 

I  loeg  +  v  T^\  "  I  /oseI  (2  Re  Toqe  7~CAj7T) 


I  (. tt')  )  ^  (2®?  *  I  ~7~C£>jp") 


(6.23) 


The  last  three  terms  in  brackets  in  (6.23)  are  the  result  of 
adding  the  (  A,  p  )  amplitude.  In  general  T(AjpN>  and  T^-rf) 
differ  in  phase  by  180  ,  and  the  term  2 Re  7~(Aj ^  T(JSj^±s  negative. 

This  is  illustrated  in  fig.  46  for  the  spectrum  of  G8 ,  and  is  in 
agreement  with  the  finding  of  Tiator,  Weber  and  Drechsel  /63 _/. 
However  when  the  other  interference  terms  of  (6.23)  are  included 
the  addition  of  (A,p) process  can  actually  enhance  the  total 
isobar  effect,  as  is  the  case  in  the  G 7  spectrum  (fig.  44)  and 
part  of  the  G8  spectrum  (fig.  45).  In  this  case  (  T0&s 

+  I  )  is  negative,  and  ToBE  is  positive  but  smaller 

_ y  ,  —  1 2. 

in  absolute  value  than  the  negative  T(Ajii)  4a,A  +  '  I  part. 

Because  of  this  cancellation  the  (A,f>)  process  causes  only  a 
slight  decrease  in  cross  section.  In  the  G10,  G11  ,  G12  spectra 
(fig.  47,  48,  49)  T^n)  interferes  constructively  with  Tose  and 
T("a  p )  interferes  destructively  with  TOBB  .  The  effect  is  that 
(  A  n  )  causes  a  large  increase  in  cross  section  while  (A ,  p  ) 
causes  a  substantial  decrease.  Thus  the  isobar  terms  can  increase 
or  decrease  the  cross  section,  depending  on  the  kinematic  situation 
When  the  (  A  tt  )  decreases  the  cross  section  the  (Ay)o  )  can  enhanc 
the  total  isobar  effect,  even  though  T^tT)  and  Tf^p')  interfere 


destructively. 
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We  find  good  agreement  with  the  experiment  except  in  the 
cases  of  the  G11  and  G12  spectra  where  the  calculation  falls 
below  the  data,  particularly  at  small  k.  A  similar  discrepancy 
between  the  G11 ,  G12  spectra  at  low  k  and  the  EED  prediction  was 
noted  by  the  experimentalists  [_  1 1_/ .  It  is  suspected  that  this 
discrepancy  is  due  to  the  rapid  variation  of  the  cross  section 
over  the  solid  angle  acceptance  of  the  photon  counters.  (The 
calculations  are  not  corrected  for  this). 

The  measured  G10  spectrum  suggests  an  increase  at  k  ^*80  MeV 
which  requires  the  inclusion  of  the  A  terms.  However  due  to 
the  crudeness  of  the  OBE  calculation  and  the  size  of  the  error 
bars  a  definite  conclusion  cannot  be  made.  Furthermore,  for  the 
measured  G7  and  G8  spectra  the  OBE  calculation  alone  provides  a 


satisfactory  fit. 


CONCLUDING  REMARKS 


We  have  used  a  simple  relativistic  model  based  on  the  one- 
boson-exchange  model  of  N-N  scattering  of  Erkelenz,  Holinde,  and 
Machleidt  /"3^_7  and  constrained  by  the  soft-photon  theorem.  The 
model  yields  ppY  cross  sections,  for  lab  energies  in  the  range 
of  158  MeV  to  730  MeV,  that  are  in  reasonable  agreement  with 
most  experiments.  The  exceptions  are  the  Orsay  measurements 
/  5^_7  and  the  G11  and  G12  spectra  of  the  UCLA  experiment  /  1 1_/  - 
they  are  data  points  where  other  calculations  show  a  similar 
discrepancy. 

We  have  considered  the  contribution  of  the  current 

diagram  (fig.  16)  to  the  pp/  cross  section.  For  the  200  MeV 
©  =  16.3°  Triumf  geometry  this  can  be  as  large  as  10%;  for 
other  experiments  it  is  negligible.  It  is  thus  another 
correction  which  must  be  considered  in  potential-model  calcul¬ 
ations  when  comparing  with  the  Triumf  data. 

The  A -excitation  diagrams,  on  the  other  hand,  are  sig¬ 
nificant  only  at  the  upper  end  of  the  photon  spectra  of  the  730 
MeV  UCLA  experiment.  The  interference  between  the  A  terms  and 
the  background  was  found  to  be  important.  In  some  cases  the  A 
excitation  lowered  the  cross  section  and  the  A  terms  with  p 
exchange  enhanced  the  total  isobar  effect. 

We  conclude  by  making  a  suggestion  regarding  future 
experiments.  It  appears  that  in  all  the  experiments  performed 
so  far  we  have  not  been  able  to  measure  significant  model 
differences.  Going  to  higher  energies  (above  200  MeV)  may 
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prove  to  be  self  defeating  as  it  will  increase  the  corrections 
and  uncertainties  in  the  potential-model  calculations.  An 
investigation  of  Fearing  0  66_7  suggests  that  asymmetries  in  pp^ 
may  be  more  sensitive  to  the  off-shell  effects.  Of  special 
interest  are  those  asymmetries  that  are  forbidden  (by  parity 
conservation)  in  elastic  scattering,  since,  in  this  case,  the 
soft-photon  theorem  determines  only  the  leading  term  in  the  k 
expansion  74  J7  *  It  may  well  turn  out  that  experiments  with 
polarized  proton  beams  are  more  suitable  for  seeing  off-shell 
effects.  The  answer  will  have  to  await  further  calculations. 

Finally,  in  view  of  the  fact  that  a  number  of  calculations 
show  a  similar  discrepancy  with  the  Orsay  data  it  would  be  useful 
to  have  the  measurements  repeated.  This  could  be  done  perhaps  by 
running  the  experiment  at  Triumf  at  the  asymmetric  angles  of  the 
Orsay  measurements. 


*The  asymmetry  along  the  direction  i  is  defined  as 
/=).  _  ( dc.  -der-V^cJV  cr_')  ,  where  is  the  pp/  cross  section 

with  the  incident  proton  spin  parallel  or  antiparallel  to  i. 
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Figure  1.  Diagrammatic 


representation  of  eq .  (1.3). 


Figure  2.  The  general  internal-radiation  process. 


1 


a 


Figure  3.  Pair-creation  current:  a  diagram  not  included 
in  a  potential-model  calculation. 
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y 


Figure  4.  Proton  momenta  in  terms  of  the  Harvard  angles. 

The  photon  momentum,  k,  (not  shown)  has  a 
negative  y  component  and  may  have  a  positive 
(0-  9r  4=  7r  )  or  negative  (  n  <  0y  <  Zt\  )  x  com¬ 
ponent,  according  to  energy-momentum  conservation. 
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Figure  6.  Definition  of  V  .  lcQ  is  the  momentum  of  the 

limiting  gamma  ray;  ot  is  chosen  such  that  c<  kQy  =  ky; 
thus  kT  =  k  -  0(  kQ  has  no  y  component.  'Yy  is  the 

angle  between  k'  and  the  z  axis. 
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Figure  8.  The  ratio  of  experimental  to  calculated 

elastic  p-p  scattering  differential  cross 
section.  Rc  =  d  <J  /d  i7cm(exp  . )  /d  <S  /d  Jlcm(calc . ) 
at  ©cm  =  6Q°,  where  d<r /d  i7cm(calc . ) 
is  the  differential  cross  section  computed  with 
the  model  of  Erkelenz,  Holinde,  and  Machleidt  [34] 
in  Born  approximation  with  the  form  factors 
modified  as  given  in  eq .  (3.3). 
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0CM  (de9> 


Figure  9.  Elastic  differential  cross  section 
computed  from  the  EHM  model  in  Born 
approximation  with  the  modified  form 
factors  (as  given  in  eq .  (3.3)). 
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Figure  11.  Elastic  (p-p)  differential  cross  sections 
from  the  four-pole  model  of  Arndt,  Bryan 
and  Macgregor  [38]  with  no  unitarization 
(solid  curves)  and  with  1  waves 
unitarized  (dashed  curves) .  In  both  cases 
the  S  wave  was  obtained  from  the  unmodified 
Born  terms. 
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Figure  12.  Elastic  (p-p)  differential  cross  sections  from  the 
four-pole  model  of  Arndt,  Bryan  and  Macgregor  [38] 
with  the  1  waves  unitarized  and  the  S  wave 
obtained  from  experiment  (solid  curve) ,  with 
L  >  1  waves  not  unitarized  and  the  S  wave  obtained 
from  experiment  (dashed  curves) ,  and  with  model 
parameters  as  modified  by  Baier ,  Kiihnelt  and 
Urban  [32]  and  no  unitarization  (dash-dotted  curves) . 
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Figure  14. 


Photon-emission  vertex. 


Figure  15.  Some  OBE  diagrams  to  third  order 

in  the  electromagnetic  interaction. 
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Figure  16.  The  UJ  radiative-decay  diagrams. 
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Figure  17.  The  tfNN  vertex  in  the  vector- 
meson-dominance  model. 
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Figure  18.  Proton-proton  bremss trahlung  differential 
cross  section  in  the  coplanar,  symmetric 
geometry.  "OBE"  is  the  result  from  the 
OBE  external-emission  model  (chapter  3) . 
"OBE±internal"  includes  the  to  radiative- 
decay  contribution  (chapter  4),  coherently 
added  with  either  relative  sign. 
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Figure  19.  See  caption  to  fig.  18.  "HJ"  is  a  Hamada-  ^ 
Johnston  potential-model  calculation.  "SPA 
is  a  soft-photon-approximation  calculation. 
The  data  points  (solid  circles)  are  pre¬ 
liminary  results  from  Trimuf  ZT  10  J7  •  The 
HJ,  SPA,  and  data  points  were  taken  from 
fig.  6  of  ref.  31. 
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Figure  20.  See  caption  to  fig.  18.  The  u>  radiative- 
decay  contribution  is  negligible  here. 
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Figure  21.  Contributions  to  the  OBE  external-emission 
model  for  Tlab  =  200  MeV,  @3  =  ©4  =  16.3 

(Triumf  geometry) .  The  curves  show  the 
coherent  addition  of  the  indicated  meson- 
exchange  contributions. 


(j-zJS/qr/)  ^p^pfyp/pp 


92 


3-2- 


1-6- 


- 7T+7}+Cf 

- 7T+7}+cr+e 


./ 


/"X. 

/  \ 

/  _ 


\ 


/ 


/*  s 


/ 


X 


/  \ 

/  \ 

/  X~X  \ 


/ 

// 


\ 


✓ 

/ 


\  \ 

\  \. 


200  MeV 
0=16.3° 


9y  (deg) 


Figure  22.  See  caption  to  fig.  21. 
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See  caption  to  fig.  21.  "Internal"  is  the 
contribution  from  the  cj  radiative-decay 
process . 


Figure  23. 
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Figure  24.  See  caption  to  fig.  21. 
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Figure  25.  See  caption  to  fig.  21. 
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Figure  26.  See  caption  to  fig.  21. 
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Figure  27.  See  caption  to  fig.  21. 
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Figure 


28.  See  caption  to  fig.  21. 
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Figure  29.  See  caption  to  fig.  18.  The  data 
points  are  from  ref.  4. 
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Figure  30.  See  caption  to  fig.  18.  Only  the  OBE 
calculation  is  shown  here.  The  data 
points  are  from  ref.  4. 
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See  caption  to  fig.  18.  The  curve 
marked  "OBE  +  #  exchange"  includes 

the  one-photon-exchange  contribution. 

The  data  are  from  ref.  48. 


Figure  31. 
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Figure  32.  See  caption  to  fig.  18. 
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Figure  33.  See  caption  to  fig.  18.  The  curve 
marked  "OBE  +  exchange"  includes 
the  one-photon-exchange  contribution. 
The  data  are  from  ref.  48. 
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Figure  34.  Differential  cross  section  for  the  Triumf 
geometry  (Tlab  =  200  MeV,  03  =  0H  =  16.3°) 
computed  from  the  OBE  external-emission  model 
The  plot  shows  the  contribution  from  the  Dirac 
coupling  (dash-dotted  curve),  anomalous 
magnetic-moment  coupling  (dashed  curve) ,  and 
the  complete  calculation  (solid  curve) . 
is  the  anomalous  moment  of  the  proton.  The 
dotted  curve  is  the  anomalous  magnetic 
moment  cross  section  divided  by  Itp  . 
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Figure  35.  See  caption  to  fig.  18.  The 
data  points  are  from  ref.  54. 
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Figure  36.  See  caption  to  fig.  18.  The 
data  points  are  from  ref.  54. 
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Figure  38.  ,A  excitation  diagram. 
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Figure  39.  ApH  vertices. 
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Figure  AO.  A  p/  vertices. 


Figure  41.  A vertex  in  the  vector- 
meson-dominance  model. 


110 


Figure  42.  Diagrams  corresponding  to  the 

terms  in  eqs .  (6.21)  and  (6.23). 
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Figure  44.  The  pptf  differential  cross  section,  d<T  /d JT^d  SZY  dk 
(nb/sr^  -  MeV) ,  as  a  function  of  the  photon  energy, 
k,  for  the  geometry  of  the  UCLA  experiment  [  11 J 

(Tlab  =  730  MeV’  @3  =  50*5°’  ^3  =  180°)-  The 

The  solid  curve  is  the  OBE  external-emission  (nucleon 
pole)  result;  the  dashed  curve  includes  the  A- 
exci tat ion-wi th- tt -exchange  diagrams  (added 
coherently);  the  dotted  curve  includes  the  total  A- 
excitation  effect,  as  described  in  chapter  6.  The 
data  points  are  from  ref.  11.  (The  angles  and 
of  this  reference  are  related  to  9*  and  'fV  , 
defined  in  section  2.1,  by  0*  =  360°  -  , 
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Figure  45.  See  caption  to  fig.  44. 
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Figure  46.  The  differential  cross  section, 
dd/dJl^dJly  dk  (nb/sr2  -  MeV) 

for  the  G8  geometry  of  the  UCLA 
experiment  (see  also  the  caption 
to  fig.  44)  from  the  A  excitation 
alone.  The  dotted  curve  includes 
only  the  7t  -exchange  diagrams  while 
the  solid  curve  includes  the  7T- 
and  the  P -exchange  diagrams. 
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Figure  47.  See  caption  to  fig.  44. 

The  long-dashed  curve  is  the 
contribution  from  the  A 
excitation  with  7T  exchange. 
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Figure  48.  See  caption  to  fig.  44. 

The  long-dashed  curve  is  the 
contribution  from  the  A 
excitation  with  7T  exchange. 
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Figure  49.  See  caption  to  fig.  44. 

The  long-dashed  curve  is  the 
contribution  from  the  A 
excitation  with  7T  exchange. 
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Figure  50.  Definition  of  the  kinematic  variables 
used  in  appendix  F. 
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Appendix  A 


The  equation  of  momentum  conservation  along  the  beam  axis 
(from  eqs .  (2.9)  and  (2.12)): 

Pi  Cos  £>3  +  pc,  COS&C,  t  (E+of  -  r3  -  E^)  cos  (9y  -  pi  ( A.  1 ) 

is  solved  here  for  =  p^  (p  ,  0y). 

Eq.  (A.1)  can  be  rewritten  as 

p<i  cos  +  P  ®  C  +  h*1'')  cos  9^  (A. 2) 

where  A  =  (  £T+0-f  ~  £~  A)  “^^Py  ~  (f  p,  -  p3  coj  @3  ; 

Squaring  (A. 2)  gives  the  quadratic  equation 

pL,  ( Co9 Qn  ~cold +  2  pq  f\  cos9if  +  (fi7--™2  cos^di)  ~  0  (A. 3) 

2  2 

When  cos  ©^  =  cos  ©y  the  solution  is 

P  =  -  (V-  (A. 4) 

2  COi  @4 

p  p 

when  cos  ©^  cos  ©*  (A. 3)  has  the  two  roots 

pty  =  -/)  cos 9h  —  [~A~  (A. 5) 

CoS^Oif  -cqs^Pv 

where  the  discriminant 

A  =  ft*  cos*  By  ~  (  cos2^v  “  co*0^)  ( A1  '  yy\roo%L6^) 

=  Coc  6y  [”  P?  +  f  to  S?@v  ”* 

The  root  with  the  lower  sign  is  in  fact  an  extraneous  root 
formed  when  (A. 2)  was  squared.  To  demonstrate  this  reexpress 
the  solutions  (A. 5)  in  terms  of  Pif  using  eq.  (A.1);  this  must 
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lead  to  an  identity  (this  procedure  is  similar  to  substituting 
the  solutions  back  into  eq.  (A.1)).  We  find 

P*  ~  L-i-i  cos  9  y  ~  cos, 

A  =  cos9y  ( ^4  UDzQi-i  ~  Ptj  c 


Pm 


(  P  Cos  Qtj  — 

\ 


Pw 


cosGy  P  )  /  ( tos*  &n  ~  cos9r 
upper  sign 


L 


'2  b-Lj  cocQc,  Cos9y  +  Pv  CcoS2Q<-i  +  CosSt) 

Cos2  &Lj  -  CO  S.Z  G  Y 


lower 

sign 


This  shows  that  the  root  with  the  lower  sign  is  not  a  solution 
of  (A.1).  The  complete  solution  of  (A.1)  is 

(P  yvP  Co  S  ■|or  OOS  Oy  ~  COS  Ocj 

2  h  Cos 

—  PccsQl y  +•  OosBy  [  /)  *  +  WpicoSpB^  ~  to  S 
co ~  to s z 9 y 


Pm 


X 


(A. 6) 


-for  cos' 9 r  ^  COS^Oij 


Appendix  B 


It  is  shown  here  that,  in  coplanar  geometry,  when 


Pi  _ J _  <<  I 

W'i  S  i  n  (  &2  + 


(B.1) 


the  kinematical  locus  in  the  plane  is  an  ellipse. 

Condition  (B.1)  implies  that  we  may  use  non-relativist ic  kine¬ 
matics  : 


P^> 


Si y\  Q 3 


pi  5 1  y\ 


^  '  k 


S  i  ^ 


O 


(B.2) 


(B.3) 


k  =  ( p,2  -  Ps  -p^VC^y^)  (BA) 

where  9g  =  ©^  +  0^,  ©^  =  9^  -  9r 

04Y  =  e4  +  ey 

Eq.  (B.4)  implies  k  <  Pi  _P|_  ^  Pi 

m 

Therefore  as  a  first  approximation  we  may  set  k  to  zero,  and  get 
an  improved  solution  by  iterating  eqs.  (B.2)  to  (B.4) 

=  0  (b.5) 


P> 


5  i  n 


0, 


(B.6) 


5 1  in 


a . 


P< 


S/  m 


0- 


(B.7) 
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Ci2> 


P< 


2  no 


(  I  -  sthlQ  i  s -77/9 


(B.8) 


si'^0. 


where 


^  "  I  5th  9*,  +  S  I'h^Pty 


5  /  h  (9,s 


Ci"') 


P3 


St  n  9< 


Pp,S(h0tj  /|  /9  S/h  Qyy') 


(B.9) 


=  pi 


Sin  9uf 
Si  to  P< 


P' 


vno 


A) 


■s  i  h  i9« 


5 1  h  P 


^  t 


cp 


p*^  p|  f  Si  n  —  Pi 


/0  5/nO  (B.  10) 


Sin  P« 


yy\ 


5 1  y\  9s  2 


where 


j/°  ^  _El 

2  no 


-  — —  an  r^in2^  +  S/Vp^v 
5<Vp0s  L 


-ft-  f\  (  5 1  h  9^  Si  Vo  P3  T  +  <S7  h  9<y  StyyOif  y'') 

yrs 

+  P1  p\2  (  Siy?  Oiy  +  Si*Z  9<ir')'^ 


HvW1 


~h 


2  no 


Pi  J_ 

wi  s>^Bs 


ft B  ~  P?  c 


m  ■ 


H 


B  “  (^/lo  Pa,  siy\93y  +  5  /  n  Pty  5/n  9^  )/  5,y ^9S 
C  =  (  Siy?  & 3Y  4  SI^Z^H O/ S/h^s 


( B .  1 1 ) 


fO 


Pi  =  P 


5/  h  9 H 
.  S>v\  9  s 


P< 


no 


-  _Pl  -1— 

^  St  v\Q$  ~3~ 


Pt  S<y>@tjY 


Pi  8  s  iy\9nr 


i  S/h*9s  2. 


Pi 


3  \  -| 


(B. 12) 
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P9 


p<  f  s/ h  e3  _  p,  _l 

L  ^  —  c/. 


[3 

Si  n  Os 


-  Pi  _i—  A  Si  VO  P 

>"V1  ^ 


3r 


~  Pi  / 


rvp  sin20s  2. 


Sih  037  +  ^  P  Pi 


Vt-A  : 


(B.13) 


('f) 


2 


p/2  F\  f  1  +  _£/  _t_  8  4-  p*  _J_- 

"*  ^  S/^9s  —  ^ 


3 

*  [B  ~  ( +■  siv^9^^)  ft  4  0(--~~z  ~77 


H 


22*  sin* 6*  )  \  (B.14) 


The  iteration  evidently  generates  a  power-series  expansion  in 


PP  S/Va  Ps  for  p^,  p^,  k*.  Therefore,  when  condition 

(B.1)  is  satisfied  the  iteration  is  convergent.  Equations 
(B.12),  (B.13),  (B.14)  give  the  first  three  terms  in  the 


expansions . 

The  kinetic  energies  T^,  T^  are: 

Ta  - 

P3 2 

"2  m 

_  p,  r  si^Otf 

o2  —  St  y~\  6^ 

p,  /  f)  Sin,  &Lj  Si  n  Oi/y 

4P  Sin  9s  SlV)6s 

(B.15) 

+  0(  1  V 

P,2  r  siv?9i 

2  m 


(B.16) 


♦The  term  .  is  included  in  the  expansion  parameter 

sin  Wg 

because  for  small  6S  high-order  terms  can  become  large,  even 

though  ^2  is  small, 
m 
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Keeping  the  first  two  terms  in  the  expansions  B.15,  B.16  the 
locus  (©■*- ) 

=  T ^  (Qy  )  is  now  shown  to  be  an  ellipse  by  expressing 

2  2 

it  in  the  standard  form  ~  =1. 

bT  q 


Let 


where  D  = 


y 


II 

T  -TC°' 

'3>  '3 

-  -D 

Si  in  6^i-y  Si  m 

(B.17) 

3t  ~ 

1 

si 

0 

=  ~D  : 

r/‘m  $3  S<  >n  ^3Y 

(B.18) 

= 

Pl’  Pi 

£ 

^  KtO  /YA 

S.Vi1^ 

-  ~D 

>S  1  m  CosOy 

5/  in  6^iy  Cos  0 y  S’  /  Jn 

(B.17) 

f  tosQy  ~  5/ h  $3  C0S&2 


(B.18) 


*  .  Si  ^  o'*  3 


0i  la  ~  5in 


'*'9^  ZC  —  D,Sf'n  dt  (  Si  va  2  $3  5/  (9<y  Cos  9^ 

+  Sl\r?&n  Siv\  9$  cos  9-$^) 

~  ~  D  siv\  6r  'Siy'&s  @i  si\n9s 


(B.19) 


Sim  0-^  Cos  ©3  cj  +  si  h  (9y  cosQi  Oc 

—  —  L)  cos&t  (  Siv\6i  Oos&i  siv\*9tf  +  Si  *9^  cos  s  i^9^)  (B#20) 

-  -  D  casQ-g-  51^03  5i'h(9ty  S/^ftr 


Therefore  from  (B.19): 

Cos^Gy  -  / 


(si*29±l 1  -  s/'tt'foy  aQ 

DZ  S/V>^@5  Si^Os 


(si*  <9,  605  Si  y  +  ft,  Cejft,  3l)  ~  c  S.lo  (9y 


and 
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coS%  ,  ^(s,se3  u,s*e3  *s;^83) 

*  ^  s'h  &3  c*>s03  S'nflv  <^1(3,  ~  s,'h*<93  s/ha£>v') 

O  S,^03  Si\*Qs 


OR 


x2  6,'SO,  +  y7*/V,203  *-  2a^s,‘hi93  s,v,0y  corftj 

—  D2  5*»  *0 3  s iv\2&n  sitr,  2 6s 


(B.21) 


define  rotated  coordinates  x',  y» 

~  *-t  ^  C-O  S(S  —  Oj  ^  S  /  h  S 

lq  ~  ocf  Si  h  S  +«j/coj<5 
and  choose  the  angle  8  so  that  the  coefficient  of  x’y*  is  zero 

.  *  .  ~  tos  8  Si  h  8  Si  \n  0Lf  +•  2  S/y£cos  8  sn^2Q 3 

+  2  (cos  8  —  siv^2  5  sty\  6  <-f  co  s  9j  —  0 

i .  e . 

Cot  2  s  ~  (sho29t/  -  Styi'lQ^)/(2  si  Vs  $l  S/vxfy  4>s.£>s)  (B.22) 

In  the  rotated  system  the  locus  is 

^  2  (  tes2<£  st'v^Bij  +-  si^S  si^2&3  +  s<*\3-S  s/yyB-j,  s '/'h<9y  cosB^) 

^  y  tt?  i  <5  Si  h  B  ^  S>  y\  S  S  t  y\  9<^  —  5/  h28s/h  By  5  /V)  Cos 

—  D2  s<'^2& 3  st'^Sif  St^lSj 


(B.23) 


Therefore,  the  ellipse  is  centred  at  T^°\  T^°^;  is  rotated  by 
angle  8  given  by  (B.22);  and  has  semi-major,  semi-minor  axes: 


a  - 


O  S i n  6A  5 1  v>  9li  si  y\  9s 


(B.24) 
'A 


(cc?52S  Sm2$iy  4  Si'h25  si*2 & 3  t  Siy\  d3  SunBy  CosB 
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L  = 


£)  S  /  v-\  0^  S/  n  0cy  Si  h  05 


(B.25) 


(  cos  S  s<  ><\  @"1,  ■*”  Jin^S  S>  *o  0^  5<  vi  2.5  5< n  S / ^  0<y  c©i 0^'')  x 


respectively.  Note  that  as  9  — ^  -ZL  (elastic  limit)  the 

S  2. 

eccentricity  goes  to  zero;  when  ©  =  I  ,  a  =  b  =  0. 

s  x 

want  a)  b,  then  (B.24)  and  (B.25)  imply  that 

7T  —  2  S  i:  0 

This  defines  the  branch  of  the  cot  function  in  (B.22)*. 
For  the  symmetric  case,  ©^  =  ©^  =  ©*  we  have 

co  s  2  0  ___ 

/)  =  2cosx&  Dsm  &3  S/V,  ~  _iL.±l_  — 

W  COS^& 


If  we 


(B.26) 


(B.27) 


T, 

p, 

cot  2 e 

(B.28) 

2  fz 

YY\ 

iS/Vn  6  (JDS  & 

Tj 

b  "  ajs 

?>_ 

rr\ 

cot  2Q 

cos *  d 

(B.29) 

b/a  " 

L&v\ 

e 

(B.30) 

i.e.  the  eccentricity  increases  as  the  opening  angle  ©  decreases; 
for  6  near  ^  the  ellipse  is  nearly  a  circle  (as  noted  before), 
and  as  9->0  (subject  to  condition  B.1)  the  ellipse  degenerates 

to  a  line  segment. 


*The  other  distinct  solution  of  (B.22) 
with  respect  to  this  solution. 


is  a  90° 


rotation 


Appendix  C 


A  Summary  of  the  QBE  Model  of  N-N  Scattering  of 

Erkelenz,  Holinde,  Machleidt  /3^+/ 

EHM  take  the  Feynman  amplitude  corresponding  to  the  exchange 
of  a  number  of  mesons  as  the  Born  approximation  to  the  full  T 
matrix.  Specifically,  if  ”p  j'p^is  the  initial,  final  CM  mom¬ 
entum  of  the  nucleons,  (  *p  [  )  is  the  relativistic  QBE 

amplitude,  and  V(~p  '  p  )  the  p-space  potential,  then  the  authors 
set 

=  (l^p)  V°8E  ^  J  ^  (~ErF) 

The  T  matrix  is  obtained  by  putting  VCp',]?  )  into  the  Lippmann- 
Schwinger  equation  (LSE) : 

TCp',  p)  =  V('p',‘p')  -  J"  ( 1  (c.2) 

Or  ) 

- 2  _  _ 5,  1  _  _J _ _ 

In  elastic  scattering  p  -  p  “  H  5 

The  LSE  is  solved  in  the  helicity  representation,  /  J  M 
then  transformed  to  the  I J  M L  S>  representat ion.  The  phase  shifts 

are  then  obtained  from  elastic  unitarity. 

The  model  parameters  are  determined  to  give  the  correct 
deutron  binding  energy  and  a  best  fit  to  the  nucleon  phase-shift 
analysis  up  to  330  MeV.  The  authors  list  the  predicted  phase 
shifts,  low  energy  and  deutron  parameters,  and  nuclear-matter 
parameters.  With  the  exception  of  the  deutron  D-state  prob¬ 
ability,  which  is  about  a  factor  of  2  too  large,  the  agreement 

with  experiment  is  good. 
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The  meson  nucleon  couplings  used  by  EHM  (in  writing  down  Vq^j?) 
are,  for  scalar,  pseudo-scalar  and  vector  mesons  respectively 

A  "  T-i  4>(s) 

3ps  T  iTs-f  <f>CeA  (C.3) 

4  =  g*  ?  TfA  T  ^  a**  f(^  4>?-  x  fi?) 

V  KV1 


For  isovector  mesons  the  fields  <^>  are  replaced  by  the  isovector 
and  the  isovector  T  is  sandwiched  between  the  nucleon  fields  ^ 
and  X  .  ^OBE  ^  ,p*  )  is  formed  by  summing  the  Lorentz-invariant 
amplitudes  of  the  various  one-boson-exchange  diagrams;  i.e., 


'N  Pi  I  Vose  P i  Pa-'? 


=  1 


<*  =  5 


f)(oih  (^3  itz') 


d  =  ps_ 


oi  -  v 


^ +  i  -pc*o  J_ _  ^ p'')  bt, 


9*0 


(C .  4) 


guo  •s'v' 


x  (-g^VtV'Ss  ■  “'"*■‘0 


. 
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where  qf  =  p4s  -  p|  =  -(p|  -  p*  ) 

and  the  subscript  on  the  spinor  denotes  both  the  spin  projection 
and  momentum. 

From  the  theory  of  integral  equations  it  is  known  that  in  order 
that  the  Born  series  from  (C.2)  converge,  the  equation  must  be  of 
the  Fredholm  type;  i.e.,  the  kernel 


K( 


P;  6 


—  Yy\ 

TSt}3 


must  satisfy  the  condition 

J  I  K^p', 


This  requirement  is  satisfied  in  the  model  through  the  use  of  form 
factors:  g^  is  replaced  by  g(oq  (qA  qM  )•  The  "coupling 

constant"  is  then  defined  as 


(C.  5) 


The  exchanged  mesons,  their  masses,  coupling  strengths  , 
goa),  and  form  factors  as  used  in  the  EHM  model  are  listed  in 


Table  C.l  below. 
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Table  D.l.  Parameters  of  the  EHM  model. 
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Appendix  D 

Phase  Shifts  from  Geometric  Uni t arizat ion 


The  partial-wave  projections,  o<  ,  are  defined  in  terms  of 
the  scattering  matrix  S  by 


—  /  4  2  C  6L 


(D.1) 


These  are  the  ex's  of  SYM  /~4o_7  divided  by  2i.  For  uncoupled 
states  the  singlet  and  triplet  amplitudes  are  and  oZjj 
(i.e.  Ujij  with  1  =  j)  respectively.  For  coupled  states 

o< 

©< 


oC  \-i  ; 


\ 


V  ^ 

-1 


the  unitarity  condition  SS  =1  reads 

ot  j  —  C  + 

Sort  ll^ll  ~  (  %rrx  H  +  ( H*  ll  ^jlO 


CD.  2) 


or 


€Xp  l  S^)  s>  n 


*Ai 

,3 


2-Xp  (l  Snry  Sjj 

(2i)~]  [  €j  exp(jc  SjLj')  -  0  7 

-i  siy\2€j  4Xp  ( i  +  <- 


=  J±I 


(D.3) 


The  (real)  Born  amplitudes  as  given  by  the  meson-pole  model  are 

denoted  by  B's,  with  the  indices  defined  as  for  the  <*»s.  We 

model  amplitude,  B,  as  the  real  part  of 


unitarize  by  taking  the 


the  unitary  amplitude  o(  ;  the  imaginary  part  of  o<  can  then  be 
obtained  either  from  (D.2)  or  (D.3).  Thus 
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8j  =  eo.s£jSiVi£j  =  i  Si'n  5j-  (D.4a) 


2,0  - 

i 

2-  <Sf 

Sji 

(  D  .  4b  ) 

J  j 

_  1 

2 

cos  2  £j 

(D.4c ) 

J 

) 

-  2. 

<^05-  6j  5 /h  *?  Sj  +  ijj 

(D.4d) 

__  1 

2 

Sir \<2  6.j  Cos  Aj 

(D.4e) 

where  Aj  ~  j  +  4h';J 

Obtaining  the  phase  shifts  for  the  coupled  states  requires  the 
solution  of  the  simultaneous  non-linear  equations  (D.4c)  to  (D.4e) 
for  7  Sj+jjJ  ?  and  .  We  used  the  Newt on-Raphson 

method  for  this:  define 


j 

*3  r 
2  s 


-f,  ~  2  CoS  2  ^  Si  Z'X-  ^J-Oj 

'h  cos, 2 2 

*  "i  ^22  CosC**^ 

For  the  zeroth  approximation  we  take  =  0,  decoupling  the 


(D.5) 


states : 


Bj-I  ;  J  *  5m 

8j +/^j  -  "i  ^ y<> 

=  0 


(D.6) 


let  (xj_,  y^,  z^)  denote  the  zero  of  (D.5)  after  i  iterations. 
By  the  Newton-Raphson  method 


where 


4b  a 
fa/ 

*C4. 


- 

Sa: 

•P, 

— 

£2 

1 _ 

1  f  'O 

X  4*, 2 

X  f\>Jj  ^3,2 

~  3c  i  +  Soc  }  4  S3  ,  2-t  +  y  -  zt  +  Sz 


(D. 


The  derivative  matrix  on  the  left-hand  side  of  (D.7)  is 


Cos  2^  cos2.o<. 

0 

~  ~2.  S/b\3-2 


o 

cos  cosZ^j 

"z.  <S/l^*?2  vS/toTa; J 


S 1  1^  5  2  s  1  ^ 

-  s/*n  <2  2  JvV\  2^ 
cosi?  2 


Appendix  E 


Quark-Model  Constraint  on  g* 

The  purpose  of  this  exercise  is  to  relate  the  sign  of  the 
ANa  vertex  to  that  of  the  NN tt  vertex  in  the  (  A,  7T )  diagram 
through  the  quark  model. 

Consider  the  NNtt  derivative  coupling: 


£. 


p-nnn  . 


ttAJM 


7T 


%  TT L  ^  ?s  If  0  'f'f 


(E.1) 


Reducing  non-relativist ically : 


~  -  jfVrfU  <T  1  IT  i.  (4, 

m7T 


By  analogy  we  take  the  non  relativistic  c^/T  coupling  to  be 


fHt  =  -J^V8  ^  S.  (E-4) 

yyi  T\ 

In  particular  consider  matrix  element  of  the  p 7T p  vertex  with 
the  pion  momentum along  the  z  direction,  and  proton  spin 
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up 

T  -  -j  O^t  Z5  cr3  i  If 3  -  .^TAjW  (E.5) 

W17I  KVItt 

We  use  the  labels  p,  n,  ?\  for  the  three  quarks  with 

<  p  I  r3 1  p>  =  l  }  <  n  I  Tj  i  o  =  -i  5  <  a  I  r3 1  ;i>  =  o 

and  make  use  of  the  quarks  structure  of  the  proton  with  spin  up: 

I  pro  lor.  -  — L  £  2  |  pT  pT  h|)?  -  /  pT  pj,  -  I  pi  pt 

+  i  I  pT  nj.  pf>  -  I  pt  »|  pi>  -  |  pj,  „1  pT>  (E.6) 

+  2)^ipTpT>  -  I  nT  pt  P!>  -/*?  pi  PT>J 

The  matrix  element  in  the  quark  model  (impulse  approximation 
with  interaction  given  by  (E.4))  is 


I  L 

yr\T{ 


pro 


(  i 


(E.7) 


Only  diagonal  terms  contribute  to  (E.7);  furthermore  the  three 
permutations  in  (E.6)  contribute  equally.  Thus 


T 


(E.8) 


Comparing  (E.5)  and  (E.8)  we  have 


(E.9) 


Now  consider  the  A  N  t\  coupling  given  in  (6.9). 
late  g*  to  ^71  $$  we  a6ain  concentrate  on  a  vertex  with 


To  re¬ 
simple 
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quark  and  spin  structure;  we  choose  the  /£+  p  p 

$2-  (  A  +)  ~  +  %  3  S2  ( pro-^n')  -  -*■  (A  . 

Using  the  identity 


A  b,;  =  /u  g_  +  ft.  b+  +  /ij  e3 


where 


A±  =  (A  it  />0//2 
B*  =  (B,  it  BO//I 


vertex  with 


(E. 10) 


with 


■n+  “  ~((f>,  ^aVj2  ~  ~  4- 

77-  -  *-4  =■ 


we  get  from  (6.9) 


£ 


AKuf  -  -  3*  Nf+f  A/  ^  TT 


where 


nc4,  =  ^Nn/u 


A 


+  ■»■  -  A0 


A//T 


/  T+-H 


4- 


A+//3  ~  A’ 


A+  4.  A°/J3 

i^+/j3  +  A‘ 


In  particular 


A++ 
ZV/TJ 


•^A”p  n1" 


=  -g*  A+t'“  3^ti+ 


The  matrix  element  of  the  vertex  in  fig.  39^.  (with  the  pion 
four  momentum  =  kM)  is  then 

T  -  ~  TZM  (  p¥)  U-C  p}  C-l 


For 


S=  = 


(pO  =  u.(V)  £**(+i) 


~Z(^i)  =  Oj  a^d  i/.ft')  J2J1> 


T 


(E. 1 1 ) 


(E. 12) 


*•?  A 


2  kv)A 


U>i  41i 
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Non-relativist ic ally 


(E.13) 


Now  we  compute  the  same  amplitude  in  the  quark  model: 


•^2 

=• 

1  Pi 

pT  pt') 

(E. 

.14) 

h  %  L+ 

v?v 

(^4-  = 

(“2  i  ■+*  4 

ib/fa) 

YV\1\ 

T  ^ 

7 

<  A++  5 

■*  -  % 

IE 

L~) 

if 

a‘rt  iV 

pro4  on 

(E 

.15) 

hon 

Only 

the  terms 

JTaf 

t 

(where 

k(~)  ~ 

=  a*rt  +  t  cr^ 

'j 

a 

contribute 

to  (E. 1 5) • 

Therefore 

T 

~  t 

1  ■**>*> 

<  A+t 

*Sz.  ~ 

if  rf  erf  Ipro-f 

OK 

m-n 

i  -  • 

(E 

.16) 

-  Hi  JT 

V^l7T 

Comparing  (E.13)  with  (E.l6)  we  get 


=  a  JI  =  _£jj 

5  m-n 


(E.17) 


VY\-;\ 
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In  the  present  case  the  two  nucleon  fields  are  on  shell 
and  the  interaction  (E.1)  can  be  rewritten  as 


uj  (n« 

o^7l  WAJ  $T\  ^  L'Ss^  ^  ^ 

3"^  ^7TN/W 

K«r>n 

(E.1 7)  and 

(E.18)  relate  g*  to  g-n  : 

9*  =  ~  5%  "4  9n  (E.19) 

With  -22-  = 

H  7T 

=  14  (E.19)  gives  _9_  —  1 . 14  x  10  5  MeV  2,  however 

47T 

we  use  the 

2 

value  of  g*  as  obtained  from  the  width  of  the  decay 

A  N  +  J\ 

=  1.875  x  10"5  MeV-2. 

H  7T 


Appendix  F 


A  Width  and  the  Value  of  g* 

In  this  appendix  the  width  corresponding  to  the  process  of 
fig.  50  is  derived  from  the  A++  p TT^  interaction  given  by 

£-  =  f  'T  W*. 

To  within  a  phase 


5  = 


—  *  —  M 


3  u-  “■  V 


J _  Vr>  I 

V3  £  ~vT  2 to. 


(?7D  8<^(jp'-p~^)  (F.1) 


the  sum  over  the  A  and  nucleon  spins  are  done  with  the  closure 


propert ies : 

^  \j 

5  p  1 


& 


P*P“  -3a--  \ 


-t-  f  Pw  P/*-')') 


and 


Y.  ~ 

5  p»>n 


=  fl 


vy\± 


(F.2) 


y~  (a  u.^  ^  ~  ~tz  ■ 


1  XA 


- 


-  j_  pt  -r,  £))<>  +  -V)  <F-3) 

>^A 


Trace  calculation: 


?*(-£%.  ?r  P*  '  >v  ? 

=  p J  p t  _L  p-p*  ”  3/*w 
A 


-  h  (p%  p*  -  pi  pA 
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■6r 

A( 

_  JL 

hoA 

(  V  pi 

=  0 

l/l 

t^A  ( 

2 

^A 

*  ¥ 
Pa 

-  3Aw  -x*-y„ 

j  ho 

(F.4) 

Pa  p*w  8  ^  0^v 


(  V  pi  -vK  pll  pj  =  -4  ( p^  pi  -  p„  pi) 


)  V"  —  — 

•  ,  Lj  Z_  LA.  LL^  LA  ^  LL 

spins 


]2  m  mA  \  kv, 


PIP  +>^ 


P/V  pp  (F.5) 

I'VU 


The  factor  ^  in  (F.5)  comes  from  averaging  over  the  initial  four 
A  spin  states. 


Decay  probabili 


ility  =  -t  I  Z  ^  ^  ^ 


^  Spins 


A 


3*  ho 

\/VW.2lo 


(2tt)V  S  (  p*  -*  p  -  $)  V T 


-2 


p-  p  +.  ™ 


VooA 


<5  hoA 


I 

I2.wi  KnA 


Replace  V"  - *  _  V  ^  £  cA  P 

^  (2  it/ 


Decay  rate 


=  p  =  it ,  9*11,  ( °*V  wws 


(2  jis^Vp*-  p-  D 


C2nr  v 


PV^W  1  ^  ^ 


X  V2  /  p-p*  vKr^ri^li)  - 

\  Wo  A  /  F  pvi 


In  the  rest  frame  of  the  A 

w=  ,  ■?=-?,  F*  = 0  and 
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6  — ! —  C  I  8  (  tnA  -  E  '  oo)  f  £T+  jvT)  2 

4n  I2tj  J  b~  cj 


Define  ZTp  =  E~  +  co 

4f  -  m2  +  ) 

=  JL  +  j_ 

60 


d  <&  =  d  Ef 

L-f 


r  =  4tiL  s^-fO^-V 

1  '  fc.P  C  60 


— 

3*  2 

Ht\ 

fm  l?l3 

3  wiA 

(F .  6) 

From 

mA  -  ( ~p  1  f  +  m/f) 

we  get 

f  +-^  ~  ^)/(2  ^ 

With 

r  « 

//V.3M>V  [b2j  )  Yn„  =  (2  32  «tV, 

Kn  =  9  /"/el/ 

and 

=  / 37  AleV 

we  get 

9*  =  /  F.7S  x  10  6  A-feV  2 

(F .  7) 

the  value  used  in  our  calculation. 


Appendix  G 

Quark-Model  Constraint  on  C 


We  write  the  non-relativist ic  amplitude  for  the  process  of 
fig.  40j  with  the  interaction  cfj  ~  -  e  Jj  C  '2TV  i  Y5 
and  in  the  quark  model.  For  simplicity  we  assume  the  A  to 
have  spin  S2  and  the  proton  to  have  spin  +-5-.  From  (6.14) 

(  Tx  )  we  have  (with  !=(+i)  =  (  \j  ij  o)  ~  “  £  +  ) 


T  =  <  C  JX  (  p~)  (  eCY)v  k^  -fev  fifxA 


tC  If  ^  k+  ~  k  £cts')  +  ) 


r$  o' 
0 


u 


(p<) 


T^P 


~  £  C  I3  Q(t  f  ^ 


(G.1) 


—  £  fj  (  ^2  ^  ( y')+V') 


Now, we  calculate  the  amplitude  from  the  quark  model. 


1  A*  Si=  +  34)  -  jf  [  I  pT  pi  nl) 

4-  |  pT  n  T  pT )  ■*"  I  h«T  pT  pT) 


(G.2) 


I  pro  T  ^  is  given  in  (E.6) . 


With  the  quark- tT  interaction 

£it>y  = 

where  ^  =  $  for  the  p  quark 

for  the  n  quark 


(G.3) 
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and  M?~  2.79  e 

2m 

The  matrix  element  corresponding  to  fig.  40^  is 


T  —  ~L/Mp  ^  T  1  €  L  0  •  (k  *  £  ( ^  A  ^2  +  ^0 

L 

—  ~L  yttp  ^  <(  pro4o*  T  |  ’Cl  CT_  (  k  *  I  A  -$2  “  + 

L 

—  2  c.  /“f5*  \ 

/y  /Ap  (  R  *  £(sV+ 

=  ~E  [  (  h  *  ^cx^\  ~~  i  ( I*  * 

“  ~  k  -X  ~  ^  ^ J  +’4  ki  ^COjjJ 

M?  t  ^COi  k+  ~  ^2  (G*4) 


2 

J3 


Comparing  (G.4)  with  (G. 1)  we  have 

£  -  -.3.  7  ?  JT 

lw 

=  -  x.\0  *  !03  Ale  V_l 


(G.5) 


This  compares  with  C  =  -2.31  x  1<T3  MeV  1  -  the  value  used  in  our 


calculation. 
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